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APRESENTAÇÃO 


A coleção de trabalhos intitulada “Coleção desafios das engenharias: Engenharia 
Civil 4” é uma obra que tem como foco principal a discussão científica por intermédio de 
diversos trabalhos que compõe seus capítulos. O volume abordará de forma categorizada 
e interdisciplinar, pesquisas cujos resultados possam auxiliar na tomada de decisão, tanto 
no campo acadêmico, quanto no profissional. 

Os trabalhos desenvolvidos foram realizados em instituições de ensino e pesquisa 
no Brasil e dois trabalhos internacionais, um localizado no México e o outro no município 
de Valenzuela, na Espanha. Em todos esses trabalhos, foram apresentadas problemáticas 
distintas que são: a inovação através do desenvolvimento de um aplicativo através da 
tecnologia Building Information Modelling (BIM) para estimar a utilização de custos de obras 
em tempo real. O uso de exemplos práticos no ensino de métodos numéricos. O estudo 
da eliminação de erros sistemáticos através de equipamentos geodésicos e topográficos 
durante a execução de locações de obras. Levantamento de causas de patologias nos 
elevados de concreto armado. Descrição através de exemplos numéricos da habilidade 
do elemento de viga unificado suportar grandes translações e rotações. E por fim, uma 
proposta de fornecer um sistema de aterro sanitário em conjunto de um plano de gestão 
ambiental que seja sustentável e economicamente viável. 

Os temas discutidos nesta obra, possuem a proposta de fundamentar o conhecimento 
de acadêmicos, mestres e todos aqueles que de alguma forma se interessam pela área 
da Engenharia Civil, através de temáticas atuais com resoluções inovadoras, descritas 
nos capítulos da coleção. Sendo assim, a divulgação científica é apresentada com grande 
importância para o desenvolvimento de toda uma nação, portanto, fica evidenciada a 
responsabilidade de transmissão dos saberes através de plataformas consolidadas e 
confiáveis, como a Atena Editora, capaz de oferecer uma maior segurança para os novos 
pesquisadores e os que já atuam nas diferentes áreas de pesquisa, exporem e divulgarem 
seus resultados. 


Armando Dias Duarte 
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CAPÍTULO 1 


DESENVOLVIMENTO DE APLICATIVO PARA 
ESTIMATIVA E ATUALIZAÇÃO DE CUSTOS DE 
OBRA EM TEMPO REAL — BIM 5D 


Data de aceite: 01/12/2021 


Dieisom Henrique Macagnan 


Engenheiro Civil, PUCRS 
Porto Alegre/RS, Brasil 


Abrahão Bernardo Rohden 

Professor PUCRS/FURB 

Porto Alegre/RS, Brasil 
http://lattes.cnpg.br/0220275442286318 


RESUMO: A aplicação da tecnologia BIM 5D na 
elaboração dos orçamentos de projeto, pode ser 
considerada como uma das últimas revoluções 
tecnológicas da indústria da construção. 
Entretendo, alguns processos ainda contam 
com a leitura e transferência de dados de forma 
manual entre diversas plataformas, sendo 
suscetíveis a falhas humanas. O objetivo deste 
trabalho foi o desenvolvimento de rotinas para 
automatizar o preenchimento, atualização e 
extração de custo dos modelos BIM. Para esta 
finalidade será utilizada a plataforma de trabalho 
Autodesk Revit 2020 em conjunto com o Dynamo 
2.1, ambos softwares da empresa Autodesk. Os 
resultados obtidos confirmam as vantagens da 
automatização do processo, garantindo maior 
agilidade e transparência durante a elaboração 
da estimativa de custo. 

PALAVRAS-CHAVE: Programação, Dynamo 
2.1, Autodesk Revit 2020, Orçamentação, Sinapi. 
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APPLICATION DEVELOPMENT 
TO ESTIMATE AND UPDATE 
CONSTRUCTION COSTS IN REAL TIME - 
BIM 5D 


ABSTRACT: The application of BIM 5D 
technology in the preparation of project 
budgets can be considered as one of the latest 
technological revolutions in the construction 
industry. Entertaining, some processes still 
rely on reading and transferring data manually 
between different platforms, being susceptible to 
human error. The objective of this work was the 
development of routines to automate the filling, 
updating and cost extraction of BIM models. 
For this purpose, the Autodesk Revit 2020 work 
platform will be used in conjunction with Dynamo 
2.1, both software from the company Autodesk. 
The results obtained confirm the advantages of 
automating the process, ensuring greater agility 
and transparency during the preparation of the 
cost estimate. 

KEYWORDS: Programming, Dynamo 2.1, 
Autodesk Revit 2020, Budgeting, Sinapi. 
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Com o aumento constante da 
competitividade entre as empresas da construção 
civil, verifica-se que o correto levantamento dos 
quantitativos e custos envolvidos no projeto 
é imprescindível na tomada de decisão de 
viabilidade ou no posicionamento perante a 
concorrência do mercado. 


Dentre as dificuldades encontradas nas 
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estimativas de custos, estão às variações entre os processos construtivos de cada empresa, 
a variação nas especificações das matérias primas e o tempo de ciclo de construção de 
um empreendimento, o que torna complexa a tarefa de gerir corretamente todas essas 
informações. Podemos, então, vincular a assertividade obtida na planilha de orçamento a 
sua qualidade, detalhamento e constante atualização da mesma (AVILA; LIBRELOTTO; 
LOPEZ, 2003). 

Para Azevedo (2016), o mercado da construção civil não tem acompanhado a 
evolução tecnológica conforme os demais setores do mercado. O autor comenta que o setor 
necessita buscar mais eficiência com o auxílio da metodologia BIM (Building Information 
Model, e da automatização dos processos em busca de mais qualidade, versatilidade e 
diferencial competitivo. 

Conforme relata Delatorre (2014), o mais novo modo de projetar, baseado na 
parametrização dos elementos, na interação entre as disciplinas durante o desenvolvimento 
do projeto e na interoperabilidade entre sistemas, ambas as melhorias de processo obtidas 
pelo uso da metodologia BIM, torna necessária a constante aprendizagem e evolução 
dos processos internos das empresas do ramo, e conceitos de trabalho envolvidos na 
concepção dos projetos. 

Em meio a essa evolução no modo de projetar, quantificar e gerenciar os 
empreendimentos, numa constante busca por uma maior eficiência e aumento de 
competitividade, evidenciamos que é necessário o desenvolvimento de processos mais 
automatizados e que permitam ao usuário um maior controle do grande volume de 
informação, visto que, no processo atual de orçamentação, ocorrem diversas entradas 
e transferências de dados entre disciplinas de projeto e plataformas diferentes de 
trabalho, permitindo que uma possível falha humana no cadastramento desses dados 
comprometa a assertividade dos projetos e orçamentos. Torna-se então, necessário um 
maior aproveitamento das informações contidas nos modelos eletrônicos utilizados para 
documentar os projetos, e a redução da transferência ou preenchimento de dados de forma 
manual entre as diversas plataformas. 

Desse modo, o objetivo deste trabalho foi o desenvolvimento de um plugin, para 
ser utilizada de modo complementar na plataforma Autodesk Revit, que permita realizar 
uma estimativa de custo direto do projeto em tempo real, baseando-se nos elementos 
constantes em modelo BIM. 

Como delimitação, não serão abordados pontos referentes a maneira de elaborar 
a modelagem, bem como não serão discutidos ou delimitados os níveis de modelagem 
adotados nos projetos por cada projetista, cabendo a eles refinar o modelo conforme o grau 
de precisão necessário. 
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2.1 BIM 


O conceito de BIM (Building Information Modeling) ou modelagem da informação da 
construção busca abranger todo o ciclo de vida de um empreendimento, sua concepção, 
projeto, detalhamento, orçamento, execução, manutenção e possíveis intervenções futuras 
de modo colaborativo e com interoperabilidade entre as diversas plataformas de trabalho e 
disciplinas envolvidas (CAMPESTRINI, 2015). 

A CBIC (2016) explica que o conceito de BIM e muito mais amplo do que aparenta, 
e pode ser visto como um conjunto de políticas, processo e tecnologias, que de forma 
integrada, configuram um novo processo de projetar uma edificação, permitindo validar o 
seu desempenho e gerenciar as suas informações em diversas plataformas com o auxílio 
de modelos virtuais. O autor afirma que esta metodologia não é nova, apenas o uso do 
termo dentro do ramo da construção civil é relativamente novo. 

Segundo CBIC (2016) devemos ter em mente que não necessariamente um modelo 
que seja realizado em um ambiente 3D pode ser considerado como BIM, para isto, tanto 
o modelo quanto à plataforma que o concebeu devem conter mais do que apenas a sua 
forma geométrica, este deve portar todas as informações necessárias sobre o ciclo de vista 
do elemento. 

Catelani (2016) salienta que dentre o mercado existem diversos softwares que 
poderiam ser descritos como BIM wash, estas são falsas plataformas de projeto que 
possuem referências 2D para simular objetos 3D, mas que não permitem a quantificação, 
atualização automática ou a realização de análises e simulações dos modelos. 

O termo BIM foi desenvolvido pelo arquiteto Phil Bernstein nos anos 90, por 
solicitação da equipe de desenvolvimento do software Autodesk Revit, classificando um 
ambiente de projeto que pudesse reunir em apenas um software diversas funcionalidades 
e permitir modelos amplamente complexos (EASTMAN ET AL., 2011). 

O mesmo autor afirma que a metodologia abrange parcelas indispensáveis do ciclo 
de vida de um empreendimento, se implementado de forma correta, agrega qualidade e 
facilidade aos projetos por contar com um sistema integrado e sem perda de informação. 
Como resultado do processo, torna-se muito mais claro verificar e solucionar interferências 
ou conflitos entre as diversas disciplinas, diluindo o processo de compatibilização, para 
que este ocorra ao longo da concepção das soluções de cada disciplina, gerando um 
produto final muito mais claro e reduzindo drasticamente as revisões de projeto por falha 


de concepção. 


2.1.1 Dimensões do BIM 


O conceito BIM possui atualmente oito dimensões: 3D, como sendo um modelo 
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virtual do empreendimento; O 4D agrega ao modelo a variável tempo, permitindo o 
planejamento de sua execução; O 5D agrega ao modelo a variável custo, permitindo a 
quantificação e levantamentos de custo; O 6D agrega informações sobre operação e ciclo 
de vida dos elementos; O 7D agrega dados de sustentabilidade para análises energéticas 
entre outras; E o 8D agrega dados sobre a segurança e prevenção de acidentes ao modelo 
(SMITH, 2014). Eastman et al. (2014) comenta que os modelos da construção possuem 
um grande potencial ainda em descoberta, e que não se limitam aos níveis já utilizados 
podendo abrangem nD's. 


2.1.2 Nível de Desenvolvimento (LoD) 


O nível de desenvolvimento trata-se de uma classificação feita pela AIA (Instituto 
Americano de Arquitetura), com o intuito de organizar as etapas do desenvolvimento de 
um empreendimento dentro do ambiente BIM. Em 2013 um documento chamado “Project 
Building Information Modeling Protoco" foi emitido pela AIA, definindo os níveis de LoD 
e especificando as suas delimitações. Os níveis são uma referência de forma clara e 
confiável sobre o conteúdo dos modelos BIM nos vários estágios do processo de construção 
conforme Figura 1 (BIMFORUM.ORG, 2019). 
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Figura 1- Níveis de Desenvolvimento (LoD) 


Fonte: http://manufacturers.bimetica.com (2019) 


Segundo CBIC (2019) os LoD são classificados conforme a Tabela 1. 


Nível de LoD Requisitos necessários 


LoD 100 Elementos representados de modo gráfico simbólico ou uma representação 
genérica. Suas dimensões e volumes são aproximados. 


Elementos representados de modo gráfico como um sistema, objeto ou conjunto 
LoD 200 genérico. São apresentadas quantidades, dimensões, formas, localização e 
orientação aproximadas. Neste nível não é necessário adicionar nenhum dado não 


gráfico aos elementos. 


Elementos representados de modo gráfico como um sistema, objeto ou conjunto 
LoD 300 específico para fins de quantidade, tamanho, forma, localização e orientação. Neste 


nível dados não gráficas podem ser adicionadas ao elemento. 


LoD 350 Adiciona-se as interfaces com outros sistemas, permitindo a compatibilização. 
Neste nível dados não gráficas podem ser adicionadas ao elemento. 
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LoD 400 Adiciona-se as informações sobre detalhamento, produção, montagem e instalação. 
Informações não gráficas podem ser adicionadas ao elemento. 


Elementos representados de modo gráfico conforme construído em termos de 
LoD 500 tamanho, forma, localização, quantidade e orientação. Informações não gráficas 
podem ser adicionadas ao elemento. 





Tabela 1 - Níveis de Desenvolvimento 


Fonte: Aluno (2019) 


2.2 Autodesk Revit 


A Autodesk (2019) classifica o software Autodesk Revit como um modelador de 
edificações em ambiente 3D com escala real, para isto baseando-se em uma estrutura 
de elementos com parâmetros específicos conforme a sua categoria ou aplicabilidade na 
edificação. 

Cada elemento que compõem o projeto possui as informações sobre sua função, 
modo de utilização, interação com os demais elementos, sua representação gráfica tanto 
em 3D como em vistas 2D e todo esse conjunto de informação é baseado em uma hierarquia 
de hospedeiros e hóspedes. Utilizando esta metodologia o modelo global da edificação 
encontra-se em constante atualização a medida que o usuário realiza uma alteração, esta 
é reproduzida em todas as representações dos objetos. A parcela de documentação do 
modelo, ou seja, suas vistas, folhas de desenho e listagens de quantitativos consistem 
apenas na apresentação das informações do modelo virtual. 

Gaspar (2015) explica que o Revit teve seu início no ano de 2000, com foco na 
disciplina de arquitetura, e Vanderzande (2013) relata que em 2013 a plataforma recebeu 
a inclusão das versões com foco em estruturas e em instalações prediais. Atualmente o 
software Autodesk Revit possui uma versão completa contendo todas as disciplinas em 
apenas uma aplicação. 

Utilizando a versão completa da plataforma, pode-se projetar e representar todas as 
disciplinas de uma edificação, desde o conceito arquitetônico e o seu lançamento refinado, 
até o zoneamento, traçado e dimensionamento das instalações elétricas, hidráulicas e de ar- 
condicionado. Toda esta informação pode ser reunida em apenas um modelo, quantificada 
e detalhada com um grande nível de precisão. 

Um importante conceito que deve ser compreendido antes de se trabalhar com esta 
plataforma são o que significam famílias, tipos e instâncias, segundo a (AUTODESK, 2019) 
as famílias são um agrupamento de elementos que compartilham sua forma geométrica 
básica e suas funcionalidades, podendo variar suas dimensões e até mesmo subpartes 
construtivas. As famílias são subdivididas em três conceitos. 

As famílias de sistema são elementos conceitualmente executados no decorrer da 
obra, sendo compostas por insumos e convertidas em um produto final diretamente no 


local. Exemplos são: telhados, forros, paredes, pisos, rampas, escadas e guarda corpo 
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(AUTODESK, 2019). 

As famílias carregáveis, ou de biblioteca, seriam elementos fabricados por outras 
empresas. Tecnicamente possuem suas dimensões e forma pré-determinada e são 
vendidos em forma de catálogos. Exemplos são: janelas, portas, mobiliário em geral 
(AUTODESK, 2019). 

As famílias criadas no local, ou “In Place” como são caracterizadas, seriam 
elementos únicos de cada empreendimento, que são confeccionados sob medida para um 
determinado espaço, e que não poderiam ser reaproveitados em outros empreendimentos. 
Exemplos são: Bancadas de suporte, mobiliário sob medida em geral (AUTODESK, 2019). 

Complementando esta estrutura, a plataforma usa o conceito de tipo para classificar 
em segundo nível os elementos, cada descrição de tipo armazena a forma geométrica, 
as dimensões e as informações que diferenciam um elemento do outro. Exemplificando, 
uma porta de madeira simples seria uma família recarregável, e estaria subdividida, por 
exemplo, em variações de tipo 80x210, 90x210, 100x210, formando assim um tipo de 
catálogo de variações disponíveis para uso no projeto. 

Por último existe o conceito de instância, onde cada elemento dentro da edificação 
pode possuir informações únicas e especiais. Um exemplo seriam duas portas com tipos 
iguais 90x210, porém com acabamentos diferentes tanto em pintura quanto ferragem, e 
situando-se em locais diferentes dentro da edificação. Deste modo cada instância de porta 


utilizada no modelo é única e pode ser classificada conforme desejado. 


2.3 Autodesk Dynamo Studio 


Segundo Autodesk (2019) o Autodesk Dynamo Studio consiste em uma plataforma 
gratuita que acompanha toda a família de produtos da empresa, podendo ser instalado 
juntamente com algum dos softwares de projeto, ou separadamente conforme o interesse 
do usuário. O aplicativo baseia-se no conceito de programação visual, transformando o 
processo convencional de elaboração de rotinas ou algoritmos de trabalho habitualmente 
feito de modo textual em um tipo de fluxograma intuitivo e visual que mantém usuário 
constantemente atualizado sobre os dados à medida que estes são processados conforme 


a Figura 2. 
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Figura 2 - Programação visual 


Fonte: https://primer.dynamobim.org (2021) 


A interação da plataforma não está somente ligada aos produtos da linha Autodesk, 
conta com interação com as linguagens básicas de programação, dentre elas, JavaScript, 
Python ou C e conexão com a maioria dos softwares atualmente utilizados. Com uma 
metodologia “Sandbox”, a plataforma permite acessar, alterar e transferir informações entre 
os softwares desejados, uma tarefa que não seria possível apenas pelo uso convencional 
dos softwares em questão. 

Em sua interface o usuário possui acesso aos chamados “Nós”, que são basicamente 
blocos de algoritmo que recebem um tipo de dado, executam uma determinada tarefa e 


fornecem dados de saída, cada “nó” é devidamente categorizado e detalhado ao usuário 
conforme a Figura 83. 
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Figura 3 - Bloco de algoritmo 


Fonte: https://primer.dynamobim.org (2021) 


Por se tratar de uma plataforma aberta e basicamente colaborativa, usuários e 
empresas de todo o mundo compartilham seus códigos e rotinas criadas mantendo um 
grande banco de dados de rotinas atualizado por meio dos portais Dynamo Dictionary, 
Dynamo Primer e Dynamo Developer Resources, entre outros fóruns e comunidade não 
oficiais. 

O funcionamento de cada “nó” é baseado no uso de algoritmos de trabalho, que 
podem ser classificados como uma sequência de tarefas e verificações que devem ser 
realizadas para solucionar um determinado problema, sendo este desde um cálculo 
matemático ou até mesmo uma ordenação ou classificação dos dados de entrada conforme 
as regras de controle previamente estabelecidas, sempre mantendo o usuário visualmente 
informado sobre o funcionamento do “nó”. 


2.4 Python 


Como complemento aos blocos de programação visual constante na biblioteca do 
Dynamo, o usuário tem disponível um nó chamado de “Python Scrip?, que permite ao 
usuário programar novas rotinas de maior complexidade e conectividade com todos os 
softwares disponíveis. 

Python.org (2019) descreve esta linguagem de programação como um processo 
interativo, interpretado e orientado a objetos, que é amplamente utilizado como extensão 
para aplicativos que necessitem de plataformas de interface programável baseados nos 
sistemas Unix, Mac e Windows. Sua estrutura é de alto nível e composta por dados e 
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digitação dinâmica organizada em classes. Esta plataforma é amplamente utilizada de 
modo amador ou comercial visto que é uma linguagem com direitos autorais, mas de uso 
aberto e com constante atualização pela organização e pela comunidade de usuários. 
Segundo Dynamo Bim (2019) a linguagem Pythone adotada por ser altamente legível, 
de fácil aprendizagem, com suporte a pacotes, módulo e com simplificada interação com 
os demais aplicativos existentes no mercado. Sua aplicação dentro do desenvolvimento da 
programação visual no Dynamo vem da simplicidade em executar avaliações condicionais 
dos dados trabalhados, permitindo substituir diversos “nós” existentes por apenas um 
código bem resumido, a Figura 4 demonstra uma rotina elaborada nesta linguagem. 
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Figura 4 - Programação em Python 


Fonte: https://primer.dynamobim.org (2021) 


31 MÉTODO DE TRABALHO 


Como ponto de partida, buscou-se identificar os processos dentro da elaboração 
dos orçamentos em plataforma 5D BIM, que demandavam mais atenção e tempo dos 
usuários. Nesse contexto, pode-se destacar como ponto crítico, a necessidade constante 
de realimentação de certos dados no software, por alteração na especificação do elemento, 


ou por uma necessidade de rever o valor atrelado a um produto ou serviço. 
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Outro ponto que foi identificado nesta avaliação, foi à necessidade de ter um 
somatório dos custos obtidos pelas diversas planilhas da plataforma, pois o software 
Autodesk Revit não permite a criação de uma planilha mestre contento todos os elementos 
constantes no modelo virtual de forma convencional, isso ocorre devida a diferença das 
categorias entre os elementos, bem como a diferença entre as unidades de medica, que 
podem ser: unidade, metro linear, metro quadrado ou metro cubico, dependendo do tipo 
de elemento. 

Tendo os objetivos traçados, foi determinada a criação de três rotinas, uma capaz 
de buscar dados em uma planilha em Excel genérica, utilizando um código mestre, sendo 
neste caso, os códigos do SINAPI, uma rotina capaz e ler as tabelas existentes no projeto, 
somar os valores dos seus respectivos totais, e apresentar ao usuário o somatório de todo 
o custo direto, e por fim uma rotina capaz de exportar os dados das tabelas para o Excel. 
Na Figura 5 foram apresentadas as principais etapas do trabalho. 
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Figura 5 - Método de trabalho 
Fonte: Autores (2021) 
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3.1 Aprendizagem da programação visual 


Na fase inicial do desenvolvimento das rotinas utilizadas no trabalho, foi necessário 
um primeiro estudo do funcionamento do software Dynamo e dos seus módulos. Para isso 
foi realizada uma série de pesquisas nos portais Dynamo Primer, Dynamo Dictionary e 
Dynamo Developer Resources, sites de suporte e de compartilhamento de informações 
entre os usuários. Dentre os compartilhamentos das comunidades, foram localizados 
alguns exemplos de aplicações de interação com o Excel e trabalho com os dados dos 
elementos do modelo, essas rotinas puderam ser estudadas e serviram de base para o 
desenvolvimento das rotinas propostas. 


3.2 Rotina para atualização de dados 


Essa rotina deve permitir ao usuário, definir um arquivo em Excel que será utilizado 
como fonte dos dados, e especificar em que coluna da planilha encontra-se cada um dos 
campos de código, descrição, unidade e custo unitário. Para isso, todos os dados devem 
ser organizados em uma lista padronizada e conter um código único para cada conjunto de 
dados. Ao executar a rotina, a mesma, busca no modelo a codificação especificada pelos 
projetistas, e com base no arquivo em Excel, alimenta os demais campos de descrição, 
unidade e preço unitário, por fim, destaca possíveis elementos que estejam sem o código 
preenchido ou por possuam um código invalido. 


3.3 Rotina para levantamento de custo direto 


Realizada a etapa de atualização dos dados, o usuário pode executar a rotina de 
levantamento de custo direto. Na tela de interface criada ele deve determinar o prefixo do 
nome das tabelas que serão quantificadas, bem como a coluna da tabela que apresenta 
o valor de custo total. Ao executar a rotina, a interface irá informar ao usuário a lista das 
tabelas que foram utilizadas no levantamento e logo abaixo será exibido o valor de custo 
direto total do modelo. 


3.4 Rotina para exportação das tabelas para o Excel 


Com o projeto devidamente ajustado, e as tabelas contemplando todos os elementos 
do modelo BIM, o usuário pode executar a rotina de exportação para o Excel. Na tela de 
interface o usuário deve determinar o prefixo do nome das tabelas que serão exportadas. 
Ao executar a rotina, a interface irá informar ao usuário a listagem das tabelas que foram 


utilizadas na exportação. A planilha em Excel criada ficará aberta na tela do usuário. 
3.5 Estudo de caso com aplicação das rotinas em um modelo 


Como estudo de aplicação para as rotinas desenvolvidas, foi tomado como 
base um empreendimento habitacional multifamiliar de quatro pavimentos com quatro 
apartamentos por andar, originalmente concebido e orçado com base em custos e serviços 
do SINAPI. Para melhor avaliação dos dados, foram utilizados dois modelos em níveis de 
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modelagem diferente, sendo um no nível de desenvolvimento LoD 200 e outro em nível de 
desenvolvimento LoD 350. Na sequência as tabelas de componentes foram comparadas 


com uma planilha de orçamentos tradicional. 


3.6 Validação com planilha de orçamento tradicional 


Para a elaboração da planilha de orçamento tradicional, foi utilizado um modelo 
proposto pela TCPO para edifício residencial de padrão médio, compatível com a 
classificação do empreendimento em questão, utilizando os custos e serviços informados 
pelo SINAPI para o estado do Rio Grande do Sul no mês de setembro de 2019. 


41 DESENVOLVIMENTO 


4.1 Pesquisa e Desenvolvimento das Rotinas 


Na fase inicial do desenvolvimento das rotinas utilizadas no trabalho, foi necessário 
um primeiro estudo do funcionamento do software Dynamo e dos seus módulos. Para isso 
foi realizada uma série de pesquisas nos portais Dynamo Primer, Dynamo Dictionary e 
Dynamo Developer Resources, sites de suporte e compartilhamento de informações entre 
os usuários. 

Durante à fase análise do processo de orçamentação, foi identificado que o software 
Microsoft Excel, é a plataforma mais amplamente utilizada na organização e trabalho com 
os dados dos empreendimentos, sendo assim, esse foi utilizado como portador dos dados 
dos elementos, e posteriormente como local de exportação dos dados das tabelas dos 
elementos do projeto. 


4.2 Configuração dos elementos no modelo 


Com o intuito de não adicionar novas variáveis ao software Autodesk Revit, foi 
utilizado o parâmetro Type Mark de cada elemento como portador do código do SINAPI. 
Para preencher o campo, o usuário deve selecionar um elemento dentro do modelo, e na 
janela de propriedades, deve editar o tipo do elemento, inserindo no campo Type Mark o 
código referente à composição ou serviço que o representa, conforme demonstrado na 
Figura 6. Esse processo deve ser executado para cada tipo de elemento constante no 
modelo. 

Se o projetista ou empresa possuir uma biblioteca de componentes, esse processo 
deve ser realizado previamente em cada um dos elementos da biblioteca, deste modo não 
será necessária à execução desta tarefa nos projetos, e será necessário executar apenas 
uma vez todo o processo de cadastro de códigos. 
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Figura 6 - Inserção do código do SINAPI em um elemento 


Fonte: Autores (2021) 


Para os materiais, foi utilizado o parâmetro Key Note como portador do código do 
SINAPI. Para preencher o campo, o usuário deve abrir o seu arquivo do projeto, ou O 
template de projeto do escritório, acessar a aba Manage, selecionar o comando Material, 
na janela de interface selecionar o material desejado na lista a esquerda, na parte direita da 
tela, selecionar a aba Identify, e por fim preencher no campo Key Note o código referente à 
composição ou serviço que o representa, conforme demonstrado na Figura 7. 
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Figura 7 - Inserção do código do SINAPI em um material 


Fonte: Autores (2021) 
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Esses parâmetros devem ser previamente preenchidos pelos projetistas de cada 
disciplina, sendo de sua responsabilidade, a correta especificação dos componentes e 


serviços utilizados no modelo. 


4.3 Rotina para atualização de dados 


A Figura 8 representa uma visão geral da rotina, cada agrupamento de cores possui 


uma função dentro da rotina. 
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Figura 8 - Rotina de atualização de dados 


Fonte: Atores (2021) 


4.3.1 Método de trabalho 


No Software Autodesk Revit 2020, deve-se abrir o projeto que se deseja trabalhar, 
e deixar ativa uma vista 3D que permita visualizar todos os elementos do modelo. Na aba 
Manage, sob o grupo Visual Programming, selecione o comando Dynamo Player, conforme 


mostra a Figura 9. 
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Figura 9 - Acionando o Dynamo Player 


Fonte: Autores (2021) 


Uma nova tela de interface será exibida, Figura 10, nessa deve ser indicada a pasta 


que contêm as rotinas salvas, e clicar em OK. 
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Figura 10 - Definindo a pasta das rotinas 


Fonte: Autores (2021) 


Neste ponto a tela de interface listará todas as rotinas disponíveis, deve-se selecionar 
o ícone abaixo da rotina de atualizar valores, conforme mostra a Figura 11. À medida que 
os campos forem devidamente preenchidos, ficarão com uma indicação de confirmação 
em verde ao lado. Estando todos configurados, pode-se acionar o botão de Run que se 
encontra ao lado do nome da rotina e aguardar a execução dela, uma mensagem de Run 


completed será exibida informando que o processo foi concluído. 
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Figura 11 - Interface de usuario - Rotina de atualização de dados 


Fonte: Autores (2021) 
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Deste ponto em diante os elementos do projeto estão com os dados atualizados, os 
objetos que foram identificados pela rotina como sem código, ou com código invalido, foram 
descritos desta maneira em seu campo de Description, permitindo ao usuário criar filtros 
de vista ou tabelas que destaquem esses elementos para posterior revisão. A rotina de 
atualização de dados pode ser executada quantas vezes forem necessárias para garantir 
que todos os elementos estejam sendo abordados e listados corretamente. 


4.3.2 Configuração das tabelas para quantificação 


Para quantificar os elementos do modelo BIM de modo correto, é necessário que 
o usuário crie as tabelas para a extração dos dados dos elementos. Devem ser criadas 
tabelas para cada uma das categorias utilizadas no projeto, bem como variações dentro de 
uma mesma categoria para cada tipo de unidade que possa existir. 

Ao executar a rotina de atualização de dados, essa identifica no final de cada uma 
das descrições a unidade que deve ser quantificado o elemento, podendo ser por área, 
volume, metro linear ou unidade. Na Figura 12 é apresentado um modelo de tabela de 
quantidades, para o levantamento de portas de modo unitário, e a Figura 13 exemplifica 
um modelo de tabela de materiais, para o levantamento de paredes por metro quadrado, 
ambas as tabelas utilizam um filtro para o parâmetro de Description buscando a unidade 
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Figura 12 - Configuração tabela de quantidade 


Fonte: Autores (2021) 
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Figura 13 - Configuração tabela de materiais 


Fonte: Autores (2021) 
4.4 Rotina para levantamento do custo direto 


Na Figura 14 é representada uma visão geral da rotina, cada agrupamento de cores 
possui uma função dentro da rotina. 
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Figura 14 - Rotina de levantamento de custo direto 


Fonte: Autores (2021) 


4.4.1 Método de trabalho 


Realizada a conferência e atualização dos dados de todos os elementos do modelo 
BIM, e a criação das tabelas necessárias para quantificar o projeto, pode-se executar a 
rotina de levantamento de custo direto, para isso, o usuário deve abrir o projeto que deseja 


trabalhar, e na aba Manage, sob o grupo Visual Programming, selecionar o comando 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 1 


Dynamo Player. 

Na lista de rotinas disponíveis, deve-se selecionar o ícone abaixo da rotina de 
levantamento de custo direto, conforme mostrado na Figura 15, e configurar os campos 
necessários até que todos fiquem com uma indicação de confirmação em verde ao lado. 
Para garantir um melhor controle, todas as tabelas que foram utilizadas nesta rotina, devem 
possuir um mesmo prefixo, garantindo que não sejam contabilizadas outras tabelas do 
projeto. Estando todos os campos devidamente configurados, pode-se acionar o botão 
de Run que se encontra ao lado do nome da rotina e aguardar a execução da mesma, 
uma mensagem de Run completed será exibida informando que o processo foi concluído 


conforme mostrado na Figura 15. 








Dynamo Player - X Dynamo Player Dynamo Player 
DN 
(2) < G < G 
Q SomaValor-RO6 SomaValor-RO6 
E 4 neady EM 4 Run completed 
AtualizaValor-RO9 | 
BA 4 ready - 
= : + Informe o Prefixo das Tabelas 4 Informe o Prefixo das Tabelas 
ExportarExcel-ROO SINAPI SINAPI 
BZ neady 


a Coluna com o Valor Total . Coluna com o Valor Total 


poi 5 
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Figura 15 - Interface de usuario - Rotina de levantamento de custo 


Fonte: Autores (2021) 


Na tela final de interface do usuário, é possível conferir a lista de tabelas que foram 
utilizadas no processo, e verificar o valor atual do custo direto do projeto. Esse processo 
pode ser realizado quantas vezes forem necessárias, para avaliar o custo do projeto, 
permitindo ao projetista realizar ajustes nos elementos em busca de um projeto mais viável 


economicamente. 


4.5 Rotina para exportação dos dados para o Excel 


A Figura 16 representa uma visão geral da rotina, cada agrupamento de cores 


possui uma função dentro da rotina. 
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Figura 16 - Módulo de exportação para o Excel 


Fonte: Autores (2021) 


4.5.1 Método de trabalho 


Com o projeto ajustado e com o seu custo verificado, pode-se executar a rotina 
de extração de dados para transportar os dados para uma planilha em Excel, para isso, 
o usuário deve abrir o projeto que deseja trabalhar, e na aba Manage, sob o grupo Visual 
Programming, selecionar o comando Dynamo Playere então selecionar a rotina de 
exportação de dados, conforme mostrado na Figura 15. Deve ser informado o prefixo das 
planilhas que se deseja exportar, acionar o botão de Run que se encontra ao lado do 
nome da rotina e aguardar a execução da mesma, uma mensagem de Run completed será 
exibida informando que o processo foi concluído conforme mostra a Figura 17. O arquivo 


em Excel criado ficará aberto na área de trabalho. 
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Figura 17 - Interface de usuario - Rotina de exportação das planilha 


Fonte: Autores (2021) 
































Coleção desafios das engenharias: Engenharia civil 4 Capítulo 1 EE 


5 | ESTUDO DE CASO 


Para a validação do devido funcionamento das rotinas desenvolvidas, foi elaborado 
um modelo BIM de um empreendimento multifamiliar, com quatro pavimentos, com quatro 
apartamentos idênticos por andar. O empreendimento em questão foi financiado e orçado 
nos moldes do programa “Minha Casa Minha Vida” com base nas composições do SINAPI. 
O modelo BIM foi subdividido em dois modelos, com diferentes níveis de desenvolvimento 
(LoD), permitindo uma avaliação do impacto da diferente complexidade dos modelos no 
custo direto total do empreendimento. 

O primeiro modelo BIM, foi elaborado com um nível de desenvolvimento LoD 200, 
contendo o modelo arquitetônico completo. Utilizando as rotinas no modelo em questão, 
pode-se obter um total de 37 itens, entre composições e serviços, correspondendo ao um 
valor de custo direto total de R$ 636.126,09. 

O segundo modelo BIM, foi elaborado com um nível de desenvolvimento LoD 350, 
recebendo a adição das disciplinas complementares de hidráulica e elétrica. Utilizando as 
rotinas nesse modelo, pode-se obter um total de 79 itens, entre composições e serviços, 
correspondendo ao um valor de custo direto total de R$ 714.983,25. 

Como verificação dos resultados, foi elaborada uma planilha típica de orçamento de 
custo direto, utilizando como referência a estrutura típica de composição de custo direto 
de um edifício comercial de padrão médio, disponível na TCPO, em sua 13º edição. O 
processo foi realizado com a extração de medidas e quantidades do modelo com nível de 
detalhamento LoD 350, e vinculação destas às composições do SINAPI em um arquivo em 
Excel. Utilizando o modelo de estrutura proposta, foi possível obter um total de 90 itens, 
entre composições e serviços, correspondendo ao um valor de custo direto total de R$ 
925.402,52. A Tabela 3 representa um resumo dos resultados obtidos. 


Fonte dos dados Itens listados Custo direto total Representatividade (%) 
Modelo LoD 200 R$ 636.126,09 





“26, 


Método Convencional R$ 925.402,52 


Tabela 3 - Resumo dos resultados 


Fonte: Autores (2021) 


Segundo a CBIC (Câmara Brasileira da Indústria da Construção), o empreendimento 
utilizado no estudo se classifica como PP 4-B (Residência multifamiliar - Prédio popular - 
padrão baixo), possui um valor de CUB (Custo Unitário Básico) de R$ 1.366,32 por metro 
quadrado de área construída, referente ao mês de setembro de 2019, para o estado do Rio 
Grande do Sul. Contando o empreendimento com um total de 924 metros quadrados, seu 
valor de custo total, segundo o índice CUB, seria de R$ 1.262.479,68. A variação de 26,7% 
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encontrada é aceitável, visto que existem diferenças entre método construtivo, tipos de 
materiais empregados e a simplificação do modelo utilizado no processo. 


61 CONSIDERAÇÕES FINAIS 


Tipicamente o levantamento do custo de um empreendimento é elaborado realizando 
medições manuais nos desenhos do projeto, transportando estes dados para uma planilha 
de orçamentos e vinculando estas medições aos índices de custo unitário. Dentre este 
processo podemos considerar que existem vários pontos de transferência de dados entre 
diferentes meios que permitem uma distorção dos dados e demandam um grande tempo 
dos projetistas. 

A automação deste processo, proposta neste trabalho, garante maior assertividade, 
agilidade na obtenção dos dados, e o controle de possíveis erros de preenchimento de 
parâmetros entre diferentes plataformas. Sua aplicação permite revisar tanto as descrições 
quanto os valores de custo de forma quase que instantânea, testando soluções e validando 
o seu impacto no custo direto do empreendimento em tempo real, direcionando o tempo 
do projetista na consolidação da melhor solução com melhor custo-benefício possível. A 
automação reduz drasticamente as atividades que não geram valor agregado, e permitem 
falhas humanas, presentes no processo tradicional de levantamento de custo direto. 

Durante a avaliação dos dados obtidos nos modelos BIM, com os dados obtidos na 
planilha convencional, pode-se verificar que o modelo com melhor nível de desenvolvimento, 
LoD 350, foicapaz de quantificar 77,3% do custo direto total obtido pelo método convencional, 
sendo esta diferença representada por 11 itens, componentes ou serviços, que não 
são tipicamente representados nos modelos BIM. Todavia, comparando os modelos em 
diferentes níveis de desenvolvimento, a inclusão das disciplinas complementares, impacta 
em apenas 8,6% no valor do custo direto total do empreendimento, representados por 42 
itens. 

Analisando os resultados, pode-se perceber que os itens faltantes possuem um 
impacto relativamente pequeno no valor do custo total direto, e que se considerarmos 
os avanços em confiabilidade dos modelos, ocasionados pela melhor representação 
dos elementos que compões o empreendimento, tais itens podem ser devidamente 
representados, resultando em uma maior assertividade nos levantamentos e ganhos de 
desempenho e confiabilidade em etapas futuras ao desenvolvimento dos modelos BIM. 

Esse trabalho reafirma que a busca constante pela melhoria no processo, 
simplificação das etapas de trabalho, e a automação de processos manuais é vital para 
nivelar a indústria da construção civil com as demais indústrias que atualmente usam 
tecnologia de ponta. 

Para trabalhos futuros, sugere-se o complemento das rotinas desenvolvidas, com 


um questionário que o projetista possa preencher com algumas informações extras, ou até 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 1 é 


mesmo extrair dados gerais sobre o projeto, buscando abranger os itens que usualmente 
não são representados nos modelos BIM, tais informações podem ser vinculadas ao projeto 


como um todo, e não propriamente aos elementos que a compõe. 
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RESUMEN: En diversos textos y cursos de 
Métodos Numéricos se presentan situaciones 
abstractas o muy artificiales alejados al área 
de conocimiento en la que se van a aplicar; es 
importante dirigir el curso en el contexto en el 
cual se va a emplear, y con esto motivar a los 
alumnos e inducirlos al mismo tiempo a las 
asignaturas específicas de Ingeniería Civil como 
Análisis Estructural, Diseho de Elementos de 
Concreto Reforzado, materias donde existen 
variadas aplicaciones de sistemas de ecuaciones 
y matrices e interpolación polinomial. En éste 
trabajo se presentan problemas de aplicación 
y se presentan conclusiones al presentar estos 
casos en clases. 
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INGENIERÍIA CIVIL 


PALABRAS CLAVE: Ingeniería, Civil, métodos 
numéricos, TecNM. 


NUMERICAL METHODS APPLIED TO 
CIVIL ENGINEERING 


ABSTRACT: In various texts and courses on 
Numerical Methods, abstract or highly artificial 
situations are presented, far removed from the 
area of knowledge in which they will be applied; 
It is important to direct the course in the context 
in which it is going to be used, and thereby 
motivate students and induce them at the same 
time to specific Civil Engineering subjects such 
as Structural Analysis, Design of Reinforced 
Concrete Elements, subjects where there are 
varied applications of systems of equations and 
matrices and polynomial interpolation. In this 
work application problems are presented and 
conclusions are presented when presenting 
these cases in class. 
KEYWORDS: Civil, 
methods, TecNM 


Engineering, numerical 


11 INTRODUCCIÓN 


Existe una preocupación en la ensefianza 
de los Métodos Numéricos en el contexto 
de la Ingeniería Civil; el profesor de manera 
general no incluye aplicaciones prácticas para 
solucionar problemas o éstos se presentan 
en forma abstracta, que para el estudiante de 
Ingeniería no tiene significado, ni mucho menos 
aplicación real, además que de manera general 
en las Instituciones de Educación Superior solo 


Capítulo 2 


existe una asignatura de Métodos Numéricos del total de los cursos de la carrera, por lo 
que es recomendable que durante la estancia de este curso se aproveche al máximo para 
conocer y dar sentido a los métodos propuestos. 

Como se describe en los estudios de ÁNGELES, L., et al (2017) del Tecnológico 
Nacional de México/Instituto Tecnológico de Altamira y en Tecnológico Nacional de México/ 
Instituto Tecnológico de Cancún por ARROYO F., CANO J., ARROYO, M. (2019), en éstas 
Instituciones se ha observado que en las carreras de Ingeniería, los índices de reprobación 
más elevados se tienen en los primeros afos, muchas veces debido a que los nuevos 
estudiantes carecen de las competencias y aptitudes para desarrollarse académicamente 
de manera óptima aunado que no encuentran sentido o aplicación real a las asignaturas; 
a pesar de obtener las puntuaciones más altas en el examen del CENEVAL (2020), para el 
ingreso del Tecnológico Nacional de México. 

Para revertir los altos índices de reprobación es ésta asignatura se realizaron 
diversas acciones para dar un significado más real a los problemas planteados en el aula y 
enfocados al âmbito de la Ingeniería Civil. Para lo cual se presentan a continuación algunos 
ejemplos desarrollados durante el curso. 

Ejemplo 1: Se presenta una situación de cálculo de esfuerzos en diferentes partes 
de una estructura tridimensional, las Fuerzas F se descomponen en 3 componentes F,, F, 


y F,; el esquema es representado en la Figura 1: 





Fig. 1 Estructura con apoyo simple para calcular la resistencia en uno de los apoyos. 


Este problema se reduce a la solución de un sistema de “m” ecuaciones lineales con 


un “n” número de incógnitas de la forma siguiente: 
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AjiXx + Aj2y +... + AjnZ = By 
As4X ER Às2y Rs AsnZ — B> 

4 + += 
AmiX + Am2y Fo + AnnZ= Bm 


Con notación matricial se escribiria de la siguiente manera: 
Aj1 Aijz AsnlyX 52] 
Avai Az Asnl|Y| | B; 
Pai And Amal'Z Bm 


O también: 
Ax =b 
Donde A es la matriz coeficiente del sistema, donde x es el vector a resolver, y b es 
el vector de términos independientes. Existen para su solución diferentes métodos directos 
como se muestran en CHAPRA, S., CANALE, R. (2015) y NIEVES,A,, DOMÍNGUEZ, F. 
(2012) como: 


* Eliminación de Gauss. 
*- Eliminación de Gauss por pivoteo. 
*- Eliminación de Jordan. 
* — Inversa por Gauss. 
*  Eilc... 
Este ejemplo se resolverá por el método Gauss-Seidel, se parte de que Ax = b para 
obtener la ecuación: 
Ax-b=0 
La ecuación vectorial correspondiente a f(x) = O, por lo que se busca una matriz B, 
y un vector C, de manera que la ecuación vectorial sea: 
x=Bx+c 
Para obtener la solución se toma, un vector inicial x? como primera aproximación al 


vector x. Posteriormente se calcula con la ecuación x = Bx + c, la sucesiones vectoriales 


x(D, x(9, x(9... x, como se muestra a continuación: 


x*+*U = Bx*+ c dondek=0,1,2....,n 
Ases efa a A o DP 


Las sucesiones x(9, x(d, x, x(9... x, deben converger al vector solución x, por 
lo que es necesario que XxX”, i<sj=<n (es decir, todos los componentes del vector x!” se 


aproximen a x, 1<j <n que corresponden a x), todas las diferencias x” - X », I<j<nsean 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 2 Bo. 


menor a un valor previamente establecido que generalmente es pequefio de acuerdo a las 
características del problema y que conserven menos los valores todos los vectores de la 
iteración: 


lim x,” = x, isj<n 


m=00 ] 

Para este ejercicio se define el algoritmo y su convergencia £; dado el sistema Ax 

= b, la manera más simple es despejar x, de la primera ecuación, x, de la segunda, y así 

sucesivamente. Se sabe que debe ser en la diagonal principal deben ser distintas de cero 
y es más rápido cuando estos valores son aproximadamente = 111. 


El sistema de ecuaciones queda de la manera siguiente: 


Aj1X + Aj2)y + Aj3W = A, 4Z — Db, 


As1X -T- As2Y -|- As 3W -[- As 4Z — b> 
Azx + Asoy + Asaw + Aggz= Ds 
Asix + Ago + Asgw + Agsz= by 


Se proponen: 


xt) =b,/A,1 ; vt =b,/A,2 , wWO=b,/A33 z(0)=b,/Ass 


Para verificar el valor e que se tiene, se sustituye en la ecuación original como se 
muestra a continuación: 


AXO + A, 27 + A, aw(9) 4 AjazO = b, O 
AXO + 4,290 + A, WO + 4,420) = bo) 
As x + As2)9) + Asa wO) + As 420) ba (O) 
E: PRE Aug Ea Aiay? + As gw (O) A Ag 4zO) ba O 


Así los valores 


Ab=by — by O 
Abi? =bo — b, O 
AbP=ba — by O 
Abg?=by — by 
Se comparan los valores de e, sí: A po SE, po EE: ul boo <Ey 
A be? < £. Se detiene el proceso. 


Para el cálculo del segundo elemento del vector x!” se sustituye x? para este caso 
resultan: 
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(1) Db) Ay (0) Aja (0) Aja (0) by — A1,27 e Aya wtO) = Ay 4280 
qgtl=—— —— Eyre — gw s— TO 
Ai Ai Aja Aja Aja 
(1) — bb Ass (1) À23 (0) Asa (0) — b, — Az XD = Az,3w'O) = Az 426) 
pla qu. yoga 
Az2 Az Az,» Az, Az,» 
(1) — by Asa, o) As a) Ass o. ba — Ag x8) — As 2 y80) — As, 2(0) 
Wit) = —— — Oy! DgtzD—D— 
Asa Aga A33 Aga, Ass 
(1) — by Ass (1) As2 (1) Ass (1) — by — Ag x) o As2y O) Es Az 3w 8) 
= Da mm “)F —  —— 4? A 5 A ii | 
Aga Ass Asa Asa Asa 


Para verificar el valor de e que se tiene, se sustituye en la ecuación original como 
se observa a continuación: 


A, x + A, 27 - A, 3w E A 42 o b, O 
AX) + As 2y ++ As 3w + As azir e b 89 
Ag 1x0) LE A3,2y 0) Ei As 3w 6 + As 4200 = b; 8) 
AXO + 4,270 + A gw + AD = b, 2 


Así los valores 


Se pretende que: A pi? > À he” ,Ã bo” >A hos A pa” >A bo, A é > 
(1) 
Ab, 
- va ad) Es 
Y nuevamente se comparan los valores de £, sí: À Db, <e, A b> E 
A b£O <EYA bo? < £. Se detiene el proceso. 
Para el cálculo del siguiente elemento del vector x? se sustituye x!” se obtienen los 


valores siguientes: 


by — Ay) — A, gw) — Ay 28) 


(2) =. 
e Aa 
po Dead A = Ag 
di | Az 
0) o Peas — Aga! — Asas” 
A3,3 
200 — Da AsaR O — Aga AsaW 
Ass 


Y se prosigue con el proceso. 
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Para facilitar la visualización de los valores de x, y, w, & z así como los de Ab, Ab,, 
Ab, & Ab,, se elabora una tabla donde se resumen estos datos como sigue: 


COCA SE 
Do [we [amy [a] um [a | 2 | as 





Tabla 1. Resumen de las iteraciones, valores y e obtenidos del procedimiento Gauss Seidel. 


En la tabla 1 se pueden observar cómo se reducen los valores de Ab,, Ab,, Ab, y Ab, 
en cada iteración; en caso que se incrementen es por algún error de cálculo que se cometa. 
El número de iteraciones para finalizar el proceso depende del valor de £ seleccionada 


desde el início, este método numérico tiene las siguientes ventajas: 


* |. Se puede obtener fácilmente una solución a groso modo desde las primeras 
iteraciones. 


* | En cada iteración se conoce el error que se comete. 
* Se considera más eficientes que los directos para sistemas muy precisos. 
* | Se puede seleccionar la precisión que se requiera desde un inicio. 


* | Emplea menos memoria que otros métodos en caso de realizar un programa 
para el cálculo. 


También tiene ciertas desventajas: 
. No se obtiene A, niel det A. 
* | Alrealizar los cálculos la convergencia está asegurada y puede ser lenta. 


* | Elsistema de convergencia puede oscilar y por lo tanto demorar más el cálculo 
si la exactitud planteada es grande. 
Ejemplo 2. Se presenta una viga en voladizo a la que se aplica una Fuerza vertical 
F en el extremo libre y tiene una deformación d, ésta última es una función de F es decir d, 
= f(F,); se miden las deformaciones al incrementar la fuerza y se obtiene la tabla 2 como 


se muestra a continuación: 
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Figura 2. Fuerza aplicada a una viga en voladizo, y su deformación. 





Tabla 2. Fuerza contra deformación. 


Se considera que la fuerza es solamente perpendicular a la viga y ésta mantiene su 
forma al aplicar la fuerza, el problema consiste en determinar si al aplicar una determinada 
F., que deformación d, tendrá, ya que el constructor establece que si ésta rebaza la 
deformación máxima d, ., puede fracturarse dicho elemento. 

Para éste problema se empleará un método conocido como extrapolación por medio 
de la aproximación polinomial, la cual requiere solucionar ecuaciones lineales. El cálculo se 
realiza fácilmente si son pocos pares de datos < 5. Si es mayor, pueden resultar complejo. 

La solución es la siguiente: la ecuación se expresa de manera general y después se 
particulariza a la viga mencionada; se parte de una función f(x) a encontrar, expresada en 


forma tabular y es de 1º grado (ecuación de una línea recta) y se escribe: 
p()=a,(4-x,) + a (x-x,) (1) 


De donde x, y x, son los argumentos de los puntos conocidos x,, f(x,), & X,, f(x,), 
y las literales a,, y a, son los coeficientes a determinar, para encontrar el valor de a,; se 


supone x = x, en la anterior ecuación y queda de la siguiente manera: 


à = PE AO) (2) 


Xi = X4 Xo = A 


Y para encontrar el valor de a,, se sustituye en valor x = x, y se obtiene como 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 2 HE 


resultado: 


p(x,) E f(x) (3) 


dl = 
" M-Xo XM-X 


Al sustituir las ecuaciones (3) y (2) en la (1) se obtiene: 





Í(Xo) fa) (x Ha) (x — X0) 
X) = MT A des ReRra 
p(x) A 1) RE 0) = Xy pm | AO são 
Para el polinomio de 2º grado (ecuación de una parábola) queda de la forma 
siguiente: 
Px) = a (x-x,) (x-x,) + à (x-x.) (x-x,)+ a,(x-x.) (X-X,) 
Donde: 


f(xo) 


E =E———————— 
(X2 — Xo)(X2 — X1) 


Para el ejemplo de este caso se empleará el método de aproximación polinomial de 
Newton donde: 


pi (x) = f(xo) + (x — x9)f [xo X1] 
FlQu)-fíxo) 


Recordando que en este método el valor de fIxo x1] = ER 
1 O 


La ecuación de 2º grado es: 


po (x) = f(xo) + (x — x9)flXo x4] + (x — xo) (x — x) FIXO, X4, X2| 


La ecuación de 3er grado es: 


pa 009) = f(xo) + (x — xo)fIxo, 2x4] + (x — xo) (x — x4)flxo, X4,X2] 
+ (x — xo) (x — x4)(x — x2)fIxo, X4,X2,X3] 


Para visualizar más fácil la tabulación de cada uno de los términos se puede generar 
la tabla 3. 
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| X 
xo | fixo) | 
| Fl] = —— —— 
| 
RR RO Flxo:x4 X:] Hu aa) dtw md 
| E Fx X2Xs] — FX XX] 
| Flx,.xs] E o - ms | 112, 22:35] RR ON = o 
[2 | Ro ER pe fara po | 
Flxz ta] = [lxa) — f Ga) | Fix x2,%3,%4] = DR 
| X3 =X 
| (xa = Flxa xa] — Fix, 2x3] 
vw | fx) AM cuia Hxs,X3 X4] =— Mão E PRE ] 5 Flxz xa. xs] — fx, x2.x3] 
Kg) = [ÃX3 E a CO 
fixa xa] = -——— ds — Xa 
| | Hg — Xs 
[| O) Flxa x4x5] = em 
| xs | Hxs) | flxs xs] = f6es) — fa) 





| ds — Xá 


Tabla 3. Tabulación general de aproximación polinomial 1º al 3er término. 


Xi RO 4º término 5º término 
xo | fixo) 
XI lix1) 
x | fix) fix Serrada x4] Fx X2.X3,X4] — fÃxosXy,X2,X3] 
qu go Ma, X3, Kg ESSES 
tas) Flxos X1,X2 XX, Xs | 
X3 X3 
Flxo XXX] — F[Xy xo, X3%4] f om pum 
Flxix2,%3Xs,%s] BR a a La La RR a a [X4, Xp, Xg, XX ] — fÃXo, Xyo Xp, Xgo Xs ] 
LÃ tdos E [WD >> >>> 
ds — dy *s — Xy 
xa | fixa) 
Xs fixs) 


Tabla 4. Tabulación general de aproximación polinomial 4º y 5º término. 


El polinomio de 4º orden: 


pa(x) = f(xo) + (x — xo)flxo, x4] + (x — xo)(x — x)fIxo, X1,x2] 
+ (x — xo) (x — x4)(x — x5)fIXo, X1,X2,X3] 
+ (x — x9)(x — x)(x — x5)(X — XD FAX» Kyo Ko, Xp X4] 


Y Ia solución es el polinomio de 5º orden: 


ps (x) = f(xo) + (x — xo)fIxo, x1] + (x — xo) (x — x9)fIxo, X4,X52] 
+ (x — xo) (x — x) (x — x2)fIX0, X14,X2,X3] 
+ (x — xo)(x — x4)(x — x2)(x — x3)fIXo, X4,X2,X3,X4] 
+ (x — xo)(x — xy) (x — xo) (x — x3)(X — X4)fÃXO, X4, X2,X3,X4, X5] 


Después de sustituir los valores de F,, F,, F,, F., F,, F, y deben obtenerse los valores 
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de la función f(x,), f(x,), f(x,), f(x,), f(x,), f(x,) Y para este caso en particular son los valores 
de las deformaciones: d,, d,, d,, d,; d,, d, - estos resultados deben ser exactos para realizar 
una extrapolación al sustituir el valor de F., y debe de ser un valor menor al de deformación 


máxima d . para evitar una posible fractura en la viga. 


max 


21 RESULTADOS 


Al explicar el método Gauss-Seidel con un enfoque inicial de una estructura con 
apoyo simple, y en el método de extrapolación por medio de la aproximación polinomial al 
presentar una viga en voladizo a la que se aplica una Fuerza vertical F versus la deformación 
d; el estudiante de Ingeniería Civil, no solo aprende el método, sino que además se interesa 
más al poder solucionar un problema en su campo de aplicación. 

Se han reducido los índices de reprobación es esta asignatura al emplear sencillos 
ejemplos enfocados a la Ingeniería Civil. 


31 RECOMENDACIONES 


Seguir aplicando ejemplos aplicados a la Ingeniería Civil y dar seguimiento para 
verificar cuanto se reducen los índices de reprobación. 


Realizar estudios comparativos a nivel nacional. 
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RESUMO: A presente pesquisa teve como 
objetivo estudar e favorecer o processo 
construtivo eliminando erros sistemáticos como 
desaprumo, locação de paredes e pilares fora 
do lugar, através da integração simultânea 
de equipamentos geodésicos e topográficos 
na execução de locações de obras. Simulou- 
se a locação de uma obra em um local plano 
e que permitiu a sua repetição pelo número 
de andares propostos (vinte andares). As 
coordenadas geodésicas destes pontos foram 
obtidas com o auxílio de equipamento GPS. 
Após processamento dos dados os mesmo 
foram convertidos ao sistema PTL, com emprego 
do programa TRANSGEOLOCAL, o qual 
realiza operações matemáticas de translação 
e de rotações de eixos em relação a um ponto 
comum aos dois sistemas (ponto origem). Para 
proceder-se a locação dos demais pavimentos, 
primeiramente 'ocupou-se, simultaneamente, 
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dois pontos aleatórios e afastados, no interior da 
obra, com o equipamento GPS. Essa operação 
permitiu a verificação da posição de cada pilar, 
em relação ao seu eixo, nos vinte andares da 
presente pesquisa. A eficácia da metodologia 
recomenda-se: o uso de equipamentos de alta 
precisão, instalação do receptor GNSS em 
tripé e base nivelante, emprego simultâneo de 
três receptores GNSS, conhecimento da teoria 
dos erros e sua aplicação para o controle dos 
enganos e erros sistemáticos. 
PALAVRAS-CHAVE: Locação de obra, GPS, 
Estação Total, TRANSGEOLOCAL, Acurácia. 


USE INTEGRATED GNSS RECEIVERS 
AND TOTAL STATION IN BUILDINGS 
LEASE 


ABSTRACT: This research aimed to study and 
promote the constructive process by eliminating 
systematic errors as plumb, lease walls and 
pillars out of place by the simultaneous integration 
of geodetic and topographical equipment in the 
execution of works of locations. Feigned to lease a 
work on a local level and that allowed its repetition 
by the number of proposed floors (twenty floors). 
The geodetic coordinates of points were obtained 
with the GPS equipment aid. After processing 
the same data were converted to PTL system 
with employment TRANSGEOLOCAL program 
which performs mathematical operations and 
translational axis of rotation relative to a point 
common to the two systems (the origin point). 
To proceed to lease the other floors, first it was 
occupied simultaneously two random points and 
apart within the work with GPS equipment. This 
operation allowed the verification of the position 
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of each pillar, in relation to its axis in twenty floors of this research. The effectiveness of the 
methodology recommended: the use of high-precision equipment, GNSS receiver installation 
on tripod and tribrach, simultaneous use of three GNSS receivers, errors theory of knowledge 
and its application to the control of mistakes and systematic errors. 

KEYWORDS: Work rental, GPS, Total Station, TRANSGEOLOCAL, Accuracy. 


11 INTRODUÇÃO 


As ciências da Geodésia e da Topografia experimentaram ao fim do último século e, 
nos primórdios deste, uma grande evolução. Evolução esta, que se faz notar, tanto no que 
se refere à teoria, como no desenvolvimento de equipamentos e programas computacionais. 

Dentre os equipamentos e técnicas novas, merecem destaque o posicionamento por 
satélites (GPS, GLONASS, Galileu, Compass/BeiDou), a medida eletrônica de ângulos e 
distâncias (estação total), o laser scanner, as imagens orbitais e suborbitais de alta resolução. 
Por sua vez, OS softwares e os hardwares voltados para os cálculos e representações 
gráficas de trabalhos geodésicos e topográficos e, que podem ser associados às obras de 
engenharia, têm facilitado e melhorado a exatidão e a representação destas tarefas. 

No entanto, no Brasil, toda essa evolução tem sido pouco adequada e utilizada 
para implantação de edificações. Ainda, a maioria das pequenas, quanto médias e grandes 
obras de edificações continua a ser implantada empregando a tradicional e ultrapassada 
técnica da execução do esquadro, com linha e trena de fita de aço, quando não de fibra, 
sobre o gabarito. E, quando da implantação de mais de um andar a locação dos pavimentos 
superiores ainda é predominante por prumada associada ao emprego de estação total. 

Algumas empresas e/ou profissionais preocupados com a obtenção de uma melhor 
qualidade de locação e, consequentemente, com a qualidade, já vem locando obras com o 
emprego de estações totais, ainda que, somente no subsolo e térreo, na maioria das vezes. 

Como é do conhecimento, os equipamentos de GNSS (siglas em inglês que indicam 
os sistemas de posicionamento por satélites artificiais) fornecem coordenadas em um 
sistema tridimensional (X, Y, Z) cuja origem é o centro de massa da Terra. Posteriormente, 
tradicionalmente, elas podem ser transformadas em coordenadas geodésicas (latitude, 
longitude e altura elipsoidal) ou ainda em Coordenadas Transversas de Mercator (UTM), 
ou suas adaptações para áreas menores (RTM e LTM) que são coordenadas cartográficas. 

No entanto, essas coordenadas não são compatíveis com o sistema de coordenadas 
no Plano Topográfico Local (PTL), pois apresentam as distorções lineares características 
deste sistema cartográfico. 

Para que os dados geodésicos, obtidos a partir de rastreio com equipamentos 
GNSS, possam ser integrados aos dados obtidos com equipamentos topográficos é 
necessário transformá-los ao PTL. BURKHOLDER (1997, 2000, 2003) Dal Forno, et. al, 
2007), .JEKELI (2006).entre outros usam da metodologia denominada por BURKHOLDER 
de 3D para efetuarem esta transformação e compatibilização. Assim procedendo, é possível 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 3 BE 


integrar dados geodésicos e topográficos para operações no PTL. Diante do acima exposto 
é possível integrar simultaneamente equipamentos geodésicos e topográficos na execução 
de locações em diversos níveis. 

Assim, este trabalho ter por objetivo de oferecer informações técnicas que possam 
contribuir para a escolha de clientes e projetistas entre a locação com equipamentos 
sofisticados e a locação convencional. Além disso, apresentar informações técnicas 
relativas à locação de gabaritos com a utilização integrada de equipamentos GNSS e 
estação total. Discutindo aspectos executivos, facilidade, agilidade, precisão e exatidão, 
listando-se suas vantagens e desvantagens. 

Expor referência que permitam a pesquisadores e projetistas optarem de uma 
maneira técnica e econômica, pela locação tradicional e a aqui descrita. 

De acordo com a NBR 13133 (ABNT, 1994, p. 3), Norma Brasileira para execução 
de Levantamento Topográfico, o levantamento topográfico é definido por: 


“Conjunto de métodos e processos que, através de medições de ângulos 
horizontais e verticais, de distâncias horizontais, verticais e inclinadas, com 
instrumental adequado à exatidão pretendida, primordialmente, implanta 
e materializa pontos de apoio no terreno, determinando suas coordenadas 
topográficas. A estes pontos se relacionam os pontos de detalhes visando à 
sua exata representação planimétrica numa escala predeterminada e à sua 
representação altimétrica por intermédio de curvas de nível, com equidistância 
também predeterminada e/ou pontos cotados.” 


Locação de uma obra é a operação inversa ao levantamento!. O sucesso da obra 
dependerá de um correto levantamento, de um projeto bem elaborado e de uma boa 
locação. 

Existem diferentes métodos de locação, os quais variam em função do tipo de 
edificação. É evidente que há diferenças em se locar um “Shopping Center” de 450 x 300 
m de área, de um edifício de vários pavimentos de 30 x 38 m de área ou uma habitação 
térrea de 8x 12 m de área. No projeto de locação a obra estará referenciada a um ponto 
conhecido e previamente definido. A partir deste ponto, passa-se a locar no solo a projeção 
da obra desenhada na planta. É comum ter-se como referência, para a locação da obra, os 
seguintes pontos: o alinhamento da rua; um poste localizado no alinhamento do passeio; 
um ponto deixado pelo topógrafo quando da realização do controle da terraplenagem; ou 


uma lateral do terreno quando este estiver corretamente localizado. 


1.1 Locação com auxílio de equipamentos topográficos 


A demarcação da obra é feita através de pontos materializados (marco, piquete, 
estacas, demarcação no meio fio) dentro ou fora do canteiro da mesma, usando-se três 
coordenadas, duas planimétricas, e uma altimétrica. (X, Y, Z), as quais são retiradas do 
1 Levantamento consiste na obtenção, em campo, das medidas de ângulos e distâncias que permitirão, em escritório, 


calcular e desenhar a superfície levantada. A locação consiste em tomarmos os dados calculados em escritório, de um 
determinado projeto de obra, e implantá-lo no terreno. 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 3 HE 


projeto e inseridas na estação total. 

A medição das distâncias, ângulos e altura é feita pela estação total, tendo apenas 
o operado de estação que se preocupa com os pontos materializados diminuindo com isto 
possíveis erros e enganos. Material empregado: Estação total; Tripé; Trena de aço; Bastão 
e Prisma. As coordenadas são retiradas do arquivo CAD ou qualquer outro programa similar 


e passadas para o terreno. 


1.2 Erros em levantamentos 


Os erros de observação são modernamente denominados “propriedades estatísticas 
das observações” GEMAEL (1994, 2004) e são devidos: 

a) à falibilidade humana proveniente de nossas imperfeições; b) à imperfeição 
do equipamento devida à fabricação, ajustes ou características básicas do mesmo; c) à 
influência das condições ambientais ou naturais que são incontroláveis e influenciam os 
instrumentos e os processos e; d) aos arredondamentos de cálculos. Nessas circunstâncias, 
o verdadeiro valor das grandezas medidas nunca é atingido, conforme as classificações 
tradicionais podem ser definidas como sistemáticos e acidentais. 

De acordo com FLORESTAL (2009, pag. 90), são feitas medidas de grandeza - 
direções, distâncias e desníveis - para que se possa fazer a representação da superfície 
da Terra. Estas medições invariavelmente são afetadas pelo acaso, por mais que os 
equipamentos utilizados sejam os melhores do mercado e por mais cuidadoso que o 
topógrafo seja em um levantamento, resultando que o valor medido não fica igual ao valor 
real. Esta diferença entre o valor real e o obtido é chamada de erro. 

Assim, os erros relacionados às medições de topografia são classificados como: 


Erros sistemáticos; b) Erros acidentais e; c) Erros grosseiros ou enganos. 


Erros sistemáticos 


Os erros sistemáticos têm as suas causas conhecidas e podem ser, portanto, 
evitados com o emprego de técnicas especiais de observação, ou eliminados posteriormente 
mediante o emprego de fórmulas teóricas como consta na obra de GEMAEL (1994,2004). 
As principais características dos erros sistemáticos são: Ocorrem sempre no mesmo 
sentido; Apresentam o mesmo valor em medições sucessivas, e; São cumulativos. 

Eles seguem, na maioria das vezes, alguma lei matemática ou física. Os erros 
sistemáticos estão ligados à imperfeição ou calibração dos equipamentos, às imperfeições 
do observador e dos métodos usados e às influências das condições ambientais. 

Para DAL'FORNO (2007), os erros quando são descobertos podem ser quantificados 
e corrigidos mediante a calibração ou testes do equipamento, e procurando-se analisar 
e compreender as interações do ambiente com as observações efetuadas. A literatura é 
pródiga em exemplos desses erros, além de procedimentos básicos de como evitá-los e 
corrigi-los. Alguns exemplos clássicos de como evitá-los podem ser referidos: A reiteração 
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e a pontaria direta e inversa nas observações angulares (erros instrumentais) e A 
colocação do aparelho a igual distância das miras no nivelamento geométrico (erro devido 


as condições ambientais). 


Erros acidentais ou aleatórios 


De acordo com DAL'FORNO (2007) os erros acidentais ou aleatórios são inevitáveis 
e ocorrem devido a causas desconhecidas e incontroláveis. Têm magnitude desconhecida, 
mas que pode ser estimada. Observações repetidas sobre a mesma grandeza, mesmo 
depois de eliminados os enganos e os erros sistemáticos, ainda podem se apresentar 
inconsistentes. Essas inconsistências “ocorrem ora num ora noutro sentido e que não 
podem ser vinculados a nenhuma causa conhecida”?. Os erros acidentais acontecem ao 
acaso e são ocasionados pelas incertezas do ser humano, imperfeições dos instrumentos 
e pelas dificuldades dos instrumentos determinarem as influências do ambiente sobre 
as medidas tomadas. Sempre se manifestam independentes de quem observa, de que 
instrumento e de que método é utilizado. 

São erros acidentais os que não podem ser previstos e de natureza desconhecida. 
Poderão apresentar valores positivos ou valores negativos. Não há como evitar esse tipo 
de erro. 

Segundo GÓÔMEZ (2003, pág. 9) “para se lidar com dados numéricos de qualquer 
tipo, é sumamente importante aprender a diferença entre uma medida e uma contagem, 
entre um engano e um erro, e entre um erro sistemático e um erro aleatório”. 

ANDRADE (1998), afirma que o conhecimento das técnicas de obtenção de valores 
confiáveis e das “medidas de qualidade” das observações é indispensável para se obter 
observações de alta precisão. Dizendo ainda que, para que o valor de uma medida possa 
ser considerado bom ou não os indicadores estatísticos de sua confiabilidade devem se 
fazer presentes. 


Erros grosseiros ou enganos 


DALFORNO (2007) afirma que os enganos são comumente denominados de 
erros grosseiros em muitos textos ou compêndios. Entretanto, eles não apresentam 
comportamento de erros, pois não seguem nem uma lei ou regra e nem estão associados 
a processos estatísticos. Possuem um paradigma totalmente diferente das inexatidões das 
medidas. Assim, conforme DAL'FORNO (2007) Apud BUCKNER (2003) a separação em 
enganos e erros pode dirimir confusões geradas por má associação do significado das 
palavras na abstração e interpretação de um problema. 

Enganos têm como causas principais as falhas humanas, associadas com a 


negligência, a desatenção, ao treinamento inadequado e a maus hábitos. Condições 


2 GEMAEL, Camil. Introdução ao ajustamento de Observações. Aplicações geodésicas. Curitiba; Editora UFPR. 1994, 
pág. 60. 
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adversas e emoções também são possíveis fontes de enganos. A principal característica 
dos enganos são deles serem relativamente grandes, quando comparados aos demais 
erros. Como exemplos de tipos de enganos, pode-se citar: A calagem displicente do 
teodolito, da estação total ou da antena de GPS em um ponto; A inversão de dígitos em 
uma leitura; A escolha errada de um pixel em uma imagem. 

De acordo com GEMAEL (1994), até em registros eletrônicos de dados, é possível a 
ocorrência de enganos. Porém, nem sempre os enganos podem ser detectados facilmente. 
Outra fonte, merecedora de atenção especial, é a tomada de decisões sem uma evidência 
ou informação suficiente. DAL'FORNO (2007) segue afirmando que o observador é que 
tem o domínio sobre os mesmos e cabe a ele o cuidado indispensável para evitá-los, ou 
de detectá-los. 


1.3 Precisão e acurácia 


Acaracterização dos erros, conforme foi realizada sucintamente acima, traz embutida 
a noção geral de precisão e acurácia (exatidão). Porém, são inevitáveis e necessárias às 
definições e a distinção entre as mesmas, principalmente no que se refere a levantamentos 
topográficos, geodésicos e/ou cartográficos, bem como, na tomada de medidas sobre 
documentos cartográficos. 

Olhando para os alvos pode-se observar que os pontos do alvo da esquerda estão 
precisos, pois apresentam uma concentração em um determinado ponto. Porém não 
acurados, pois se encontram afastados do centro do alvo (verdadeiro valor ou verdade). 
Mas, se observar o alvo da direita verá que os pontos encontram-se acurados, ou seja, 
precisos e concentrados no centro do alvo, onde é a sua referência. 

O que se pode concluir é que precisão é a concentração de pontos em torno 
da média, e acurácia é a concentração de pontos precisos em torno de um referencial 


verdadeiro. 


2 | MATERIAIS E MÉTODOS 


Este capítulo mostra como a pesquisa foi realizada, os equipamentos empregados 
e suas etapas dentro do escopo metodológico a fim de apontar como as conclusões foram 
feitas da forma que se propõe. 


2.1 Materiais 


Para a execução deste trabalho foram empregados os seguintes equipamentos: 
Estação total e acessórios: 

Estação Total marca RUIDE modelo RTS — 820 Rº apresentada na Figura 5, com 
precisão angular de 2”, conforme as especificações da NBR 13.133/94, Tripé para apoio 
da mesma, bastão e prisma para o auxílio da visada ré direta e inversa, mini prisma para 
locação e coleta dos pontos. Além disso, foi utilizado também um receptor GNSS GPS, 
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marca SOUTH modelo S86-5. 

Para a transformação das coordenadas geodésicas em coordenadas no plano 
topográfico utilizou se do programa TRANSGEOLOCAL V2.1, disponível em hittp:// 
sevenengenharia.eng.br/downloads/ 

Para a representação dos resultados gráficos valeu-se do programa Auto CAD Civil 
3D. 

O processamento das observações GPS foram executadas no programa SSO, 
TopCon Tools. Apesar dos dados terem sidos coletados em duas frequências, pela 
incompatibilidade entre o equipamento e o software, só foi possível obter as coordenadas 
a partir da frequência L1. Os cálculos matemáticos e estatísticos foram executados no 
programa Excel. 


2.2 Local do experimento 


A presente pesquisa foi efetuada em uma área de estacionamento não ativado ao 
trafego, localizado aos fundos do posto Petrobras, dentro do Campus da Universidade 
Federal de Santa Maria. 

Como o objetivo do presente trabalho é viabilizar a integração simultânea de 
equipamentos geodésicos e topográficos na execução de locações de obra em diversos 
níveis, e como não havia disponibilidade de edificação para a execução real do trabalho, 
simulou-se a locação de uma obra em um local plano e que permitisse a sua repetição pelo 
número de andares propostos. 

Optou-se por simular uma construção constante de um pavimento térreo e vinte 
andares. Tal escolha se deve a facilidade para se estabelecer percentagem estatística 
(95% corresponde a 19 andares). 

Para o desenvolvimento do trabalho, realizaram-se a implantação de três pontos 
de base, externos a área do prédio, estáveis e permanentes, que permitissem, a partir 
deles, a locação em uma área de 420,54 m? dos eixos dos pilares que corresponderiam ao 
pavimento térreo. 

As coordenadas geodésicas destes pontos foram obtidas com o auxílio de 
equipamento GPS, com tempo de rastreio de trinta minutos, em cada ponto. Após 
processamento dos dados os mesmo foram convertidos ao sistema PTL, com emprego do 
programa TRANSGEOLOCAL. 

De posse das coordenadas no PTL, e do levantamento topográfico prévio do 
local, usando o programa Auto CAD Civil 3D, determinou- se as coordenadas dos pilares, 
georreferenciadas aos pontos de base. 

Estes pontos dos pilares, após a sua implantação serviram como testemunhas para 
a locação dos vintes pavimentos superiores. 

Para proceder-se a locação dos demais pavimentos, primeiramente ocupou-se, 
simultaneamente, dois pontos aleatórios e afastados, no interior da obra, com o equipamento 
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GPS. O tempo de rastreio em cada ponto foi de 30 minutos. Esta operação foi repetida vinte 
vezes, assumindo-se para cada operação dois pontos diferentes. 

No levantamento geodésico foram determinadas as coordenadas geodésicas de 
quarenta pontos (PV1iA, PV1B, PV2A, PV2B, PVSA, PV3B, PV4A, PV4B, PV5A, PV5B, 
PV6A, PV6B, ..., PV20A, PV20B), implantados na área estudada com auxilio de pregos 
de aço, fixados no chão com marreta e demarcados com auxílio de tinta. As coordenadas 
geodésicas, de cada par de observações, foram transformadas para o mesmo PTL do 
andar térreo, (através do programa TRANSGEOLOCAL), e a partir daí associadas às 
coordenadas dos pilares. Essa operação permitiu a verificação da posição de cada pilar, 
em relação ao seu eixo, nos vinte andares da presente pesquisa. 

As diferenças de coordenadas dos eixos dos pilares obtidas em relação ao andar 
térreo (testemunha), adquiridas em cada repetição, foram tabuladas e posteriormente 
calculadas as distâncias correspondentes. De posse dessas informações calculou-se as 
médias e os desvios padrões para o total de trinta e seis pilares de cada pavimento e, para 


cada pilar dos vinte pavimentos. 


31 RESULTADOS 


Este capítulo tem como objetivo apresentar a análise e a discussão dos resultados 
obtidos. Assim sendo, foram analisadas as imprecisões das coordenadas de cada pilar, em 
todos os andares, quando comparadas com as do pavimento térreo. 

Os dados dos 36 pilares para cada um dos pavimentos, com seus respectivos 
deslocamentos nos eixos x e y e seu deslocamento total, todos expressos em centímetros. 
Nas duas últimas linhas das Quadro são apresentadas as médias dos deslocamentos totais 
e o seus desvios padrões entre todos os pilares. 

Se verificou discrepâncias acima das aceitáveis para o objetivo deste trabalho. 
Estes resultados desprezados tem a característica de apresentarem um erro sistemático 
(concentração dos desvios no primeiro quadrante e com valores acima do dobro dos obtidos 
nos demais pavimentos). Na realidade, estes erros sistemáticos tem origem em um engano 
(erro grosseiro), ou seja, a falta de verticalidade do bastão suporte do GPS. Chegou-se a tal 
conclusão, a partir da constatação que nestes pavimentos os dados coletados por rastreio 
GPS, foram efetuados por um ajudante do autor. Tal ajudante, com formação de segundo 
grau, tem pouca informação sobre controle de erros. 

Observou-se que nestes pavimentos as médias, na maioria dos casos, são 
superiores ao dobro das médias obtidas pelo autor e, os desvios padrões também são 
mais elevados, o que caracteriza uma maior dispersão dos dados e consequentemente 
uma menor precisão. 

Para efeito de estabelecer critério de boa locação, o autor estabeleceu que seria 


recomendável um deslocamento de um centímetro no centro do pilar, e aceitável até uma 
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vez e meia este deslocamento. 

Constatou-se que no pilar P06 as discrepâncias em 95% de suas locações apresentam 
valores superiores ao estabelecido pelo autor como recomendável, chegando a atingir o 
valor máximo de 4,9 cm, no 2º pavimento. Esta elevação constante nas discrepâncias leva 
a crer que houve um deslocamento no pino de aço testemunha do pilar, não superior a um 
centímetro. 

Este fato não ocorreria em uma locação real de uma obra, uma vez que, cada 
andar seria marcado independentemente da materialização do ponto do andar térreo. 
Provavelmente este ponto foi deslocado pelo tráfego de caminhões pesados no local do 
experimento. 

Observando-se a Quadro 1 é possível verificar que as discrepâncias das 720 locações 
(36 x 20) em 64,72 % dos casos é inferior ao valor estabelecido como recomendável e, em 
75, 69 % dos casos é inferior ao valor aceitável. 


Ti ro Qro 


Quadro 1 — Distribuição das discrepâncias totais em número e percentagem. 





Fonte: O autor. 


Quando se desprezam os seis primeiros pavimentos, cujos resultados das 
discrepâncias foram considerados contaminados por erros grosseiros (enganos), as 
distribuições das discrepâncias apresentadas no Quadro1 passam a ter 98,61 % de 
resultados iguais ou inferiores ao valor aceitável (497 pontos de um total de 504). 

Para melhor visualização dos resultados dos Quadro foram feitas a representação 
gráfica da distribuição das discrepâncias nos pilares de canto. 

Observando a Figura 1, a qual representa o pilar 1 e tomando como base o pilar do 
pavimento térreo, pode-se observar que 60% dos pontos encontram-se acurados. Em 12 
andares as locações do pilar 1 estão com deslocamento menor ou igual a 1 cm e 75% está 
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com seu deslocamento menor ou igual a 1,3 cm, valor este menor do que o aceitável. 


Tendo em vista que, aplicando-se este método de locação os erros não se propagam, 


pois cada andar é locado independentemente do outro. No nosso caso o deslocamento do 


centro do pilar PQ01, no 20º pavimento, em relação ao térreo, ficou em 0,5 cm. 


PILAR 1 
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adam 


| SEM ESCALA 





Figura 1 - Discrepancias de locação, por andar no pilar PO1. 
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Na Figura 2, que representa o Pilar 3, pode-se observar o que segue: i) quando 


se consideram todos os andares as locações apresentam erros grosseiros (andares 2, 


3,4 e 6); ii) erros sistemáticos em todos os outros andares (viés para o 1º quadrante); iii) 


baixa precisão (grande dispersão) e iv) em função das constatações anteriores pode-se 


concluir neste caso houve erros grosseiros na instalação do GPS para o primeiro caso e, 


possivelmente, pequeno deslocamento do pino testemunha. 
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Figura 2 - Discrepancias de locação, por andar no pilar POS. 


Na Figura 3 que representa o pilar P34 os pontos encontram se acurados e 70% 
deles encontram se dentro do raio de 1 cm e, dos seis pontos que se encontram externos a 
este raio, cinco deles são pontos suspeitos de possuírem um erro grosseiro. 

Por sua vez a Figura 3 que representa o pilar P36, apresenta as mesmas 
características e percentagem Figura 4. O que só vem a reforçar que o operador do 
aparelho Receptor GPS, deve ter extremo cuidado, para que o mesmo fique no prumo, 
para melhor acurácia do seu trabalho. 
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Figura 4 — Discrepancias de locação, por andar no pilar P36. 


Coleção desafios das engenharias: Engenharia civil 4 


Capítulo 3 


IVOVOOVVVOVOOVOCCOVOVCOO 


LEGENDA 


m 
Es) 
as) 
m 
O 


ar 8 


o [6] 
. 


ed 


LEGENDA 


TERREO 


> e) q > o) La) = 


"o 


* ANDAR 
* ANDAR 
* ANDAR 


* ANDAR 


ANDAR 


* ANDAR 


ANDAR 


ANDAR 


* ANDAR 


* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 
* ANDAR 


* ANDAR 





41 CONCLUSÕES E RECOMENDAÇÕES 


O método proposto se mostrou eficaz para os propósitos desta pesquisa e pode 
ser recomendado para a locação de edifícios convencionais de vários pavimentos. Os 
resultados obtidos, com a sua utilização, são superiores aqueles resultantes da aplicação 
do método tradicional. 

A sua exatidão, quando se desprezam os pontos afetados pelo erro grosseiro 
(engano) fica, em 83,93 % dos casos menor ou igual ao recomendado (1 cm) e 98,61 % 
inferior a 1,5 cm, valor este perfeitamente aceitável em locações convencionais. 

Os resultados obtidos podiam ser menos discrepantes se tivessem sido empregados 
equipamentos mais sofisticados, atualizados e calibrados. A ação de dois operadores, com 
formação diferenciada, também influiu nos resultados. 

Tomando-se os cuidados acima recomendados, pode-se dizer que o método, aqui 
proposto, é adequado, acurado, rápido e não propaga erros que por ventura ocorram na 
locação do pilar em um pavimento. 

Para que se melhore a eficácia da metodologia aqui sugerida, recomenda-se: 

O uso de equipamentos de alta precisão: estação total com precisão angular 
superior a 2” e linear a 2 mm/km +2 ppm, equipamentos GNSS de duas frequências (L1/ 
L2); Instalação do receptor GNSS em tripé e base nivelante; Emprego simultâneo de três 
receptores GNSS, o que permite o emprego de triangulação e consequentemente maior 
precisão nas coordenadas dos pontos de controle; Conhecimento da teoria dos erros e sua 
aplicação para o controle dos enganos e erros sistemáticos é recomendável e de extrema 
importância; Para controle das influencias ambientais recomenda-se que a locação de 
todos os andares, seja feita, de preferencia, em um mesmo horário. Sugere-se que os 
trabalhos, tanto de coleta de dados, quanto de locação, sejam executados numa mesma 
faixa de temperatura; O emprego de funcionários altamente treinados. 
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RESUMO: Estruturas de concreto podem 
apresentar patologias de causas diversas, 
podendo ter origem na fase de projeto, na fase 
de execução das obras e até na fase de uso. 
Neste trabalho, realiza-se um levantamento nos 
elevados da Região Metropolitana de Belém 
a fim de identificar as possíveis causas das 
patologias apresentadas nestas obras. Expostas 
essas causas, são exemplificados os casos mais 
usuais, apresentando as possíveis soluções 
para os problemas encontrados. A durabilidade 
das estruturas de concreto é determinada pela 
qualidade e composição dos materiais utilizados, 
dosagem ou traço do concreto, compactação ou 
adensamento efetivo na estrutura, processo de 
cura efetiva e cobrimento nominal das armaduras. 
Os resultados mostram que, nas obras dos 
elevados visitados, problemas de manutenção 
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destas estruturas e até redução na seção da 
armadura exposta, ocasionado por processo de 
corrosão, caracterizando diminuição da seção da 
armadura. Sabe-se que a perda de seção das 
armaduras contribui, significativamente, para a 
perda da capacidade de carga das estruturas 
de concreto, podendo induzilas ao colapso. 
Além disso, foram identificadas similaridades 
nas principais causas das falhas construtivas 
encontradas na entrega pós-obra. 

PALAVRAS-CHAVE: Patologias; estruturas de 
concreto; corrosão; manutenção. 


CAUSES OF PATHOLOGIES IN HIGH- 
WORKS OF THE METROPOLITAN 
REGION OF BELEM 


ABSTRACT: Concrete structures may present 
pathologies of different causes, which may 
originate in the design phase, in the execution 
phase of the works and even in use. In this work, 
a Survey is carried out in the highlands of the 
Metropolitan Region of Belém in order to identify 
the possible causes of pathologies presented 
in these structures. Having exposed these 
causes, the most usual cases are presented, 
presenting the possible solutions to the problems 
encountered. The durability of concrete structures 
is determined by the quality and composition 
of the materials used, dosage or trace of the 
concrete, compaction or effective densification 
in the structure, effective curing process and 
nominal cover of the reinforcement. The results 
show that, in the works of the high visited, 
problems of maintenance of these structures 
and even reduction in the section of exposed 
reinforcement, caused by corrosion process, 
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characterizing a decrease of the section of the armature. It is known that the section loss of 
the reinforcement contributes, significantly, to the loss of the load capacity of the concrete 
structures, being able to induce them to collapse. In addition, similarities were identified in the 
main causes of constructive failures found in post-work delivery. 

KEYWORDS: Pathologies; concrete structures; corrosion; maintenance. 


11 INTRODUÇÃO 


O Brasil é um país em pleno desenvolvimento econômico e, mesmo assim, convive 
com diversos problemas sociais (desemprego, saúde, educação, habitação, entre outros), 
de tal modo, não pode se dar ao luxo de custear despesas extremamente elevadas em 
obras de reparo estrutural, as quais, na maioria das vezes, poderiam ser evitadas. 

O estudo de manifestações patológicas em estruturas de concreto armado, conhecida 
como patologias das construções, é um novo campo da engenharia das construções que 
se propõe a pesquisar as origens, formas de manifestação, consequências e mecanismos 
de ocorrência das falhas e dos sistemas de degradação das estruturas (SOUZA e RIPPER, 
1998). 

Os fenômenos patológicos habitualmente apresentam manifestação externa 
característica, a partir da qual se pode deduzir a natureza, a origem e os mecanismos dos 
fenômenos envolvidos (HELENE, 1992). Dessa forma, entende-se que diversas são as 
causas que levam uma estrutura de concreto armado à sofrer danos, por isso é de extrema 
importância o estudo desse ramo da engenharia, para que sejam evitadas manifestações 
patológicas que venham a diminuir a durabilidade das estruturas, assim como é necessário 
um conhecimento de como solucionar e recuperar aquelas que apresentam problemas 
patológicos, de maneira a recuperá-las, curá-las e impedir que o agente causador volte à 
causar os mesmos problemas. 

Determinadas manifestações incidem com mais constância, devido à necessidade 
de cuidados que frequentemente são ignorados, seja no projeto, na execução ou até 
mesmo na utilização. Segundo Silva (2011), pode-se dizer que os problemas patológicos 
de maior gravidade nas estruturas em concreto armado, especialmente pelo seu evidente 
risco à integridade da estrutura, são a corrosão da armadura do concreto, as fissuras e as 
flechas excessivas das peças estruturais. 

Por outro lado, a ABNT NBR 6118:2014 destaca como mecanismos preponderantes 
de deterioração do concreto a lixiviação por ação de águas puras, carbônicas agressivas 
ou ácidas, que dissolvem e transportam os compostos hidratados de pasta de cimento; 
a expansão por ação de águas e solos que estejam contaminados com sulfatos, dando 
origem a reações expansivas e deletérias; a expansão por ação das reações entre os 
àlcalis do cimento e certos agregados reativos; e as reações químico-físicas superficiais 
de certos agregados, decorrentes de transformações de produtos ferruginosos presentes 


na sua constituição mineralógica. A mesma norma, também destaca, os mecanismos 
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preponderantes de deterioração relativos à armadura, como o processo de carbonatação 
nas estruturas de concreto armado e a despassivação por elevado teor de cloreto. 

A soma de tantos fatores prejudiciais à estrutura, reforçada pelo crescimento 
acelerado que a construção civil vem atingindo, provoca a necessidade de inovações que 
trazem em si a aceitação subentendida de maiores riscos, leva ao ponto de viver em uma 
época de grandes preocupações. 

Na construção civil a rapidez na execução e o emprego de novas tecnologias são 
fundamentais para tornar o mercado mais competitivo e possibilitar novos empreendimentos. 
Inovações e rapidez de execução acarretaram na necessidade de se desenvolver o 
conhecimento relacionado às estruturas e aos materiais, a fim de proporcionar o avanço 
das construções (MARTINS E FIORITI, 2016). 

Entretanto, ainda que o investimento em tecnologias relacionadas a essas 
construções cresceu e evoluiu desde o início das civilizações até os dias atuais, ainda assim, 
existem sérias limitações e inevitáveis falhas que acabam por prejudicar o desempenho de 
algumas estruturas. Esses fatores conduzem ao surgimento de manifestações patológicas, 
que podem ser entendidas como perda ou queda de desempenho de um conjunto ou 
componente estrutural. 

Como a patologia pode ser influenciada pelo comportamento da estrutura em uso, 
pelo tempo e pelas condições de exposição, percebe-se, portanto, que há uma forte relação 
entre a patologia e desempenho, vida útil e durabilidade da edificação (SANTOS, 2012). 

A rapidez na execução e o emprego de novas tecnologias são fundamentais para 
tornar o mercado mais competitivo e possibilitar novos empreendimentos. Inovações 
e rapidez de execução acarretaram na necessidade de se desenvolver o conhecimento 
relacionado às estruturas e aos materiais, a fim de proporcionar o avanço das construções. 

Portanto, este trabalho trata de um levantamento de causas de patologias nos 
elevados da Região Metropolitana de Belém, onde através de inspeção visual foram 
identificadas as manifestações patológicas existentes nos sistemas estruturais em concreto 
armado dessas obras. Salienta-se, contudo, que não fez parte do objetivo, entrar no mérito 
da qualificação e da atuação dos profissionais, assim como das empresas que participaram 
dos projetos e execução destes empreendimentos, sendo o único foco a identificação das 
manifestações patológicas vistas sob a ótica da sintomatologia. 

Justifica-se este trabalho na necessidade de se evitar a ocorrência de manifestações 
patológicas no concreto armado, que acabam por diminuir a vida útil e desempenho das 
edificações, além de gerar um acréscimo de custo para recuperar e reforçar uma estrutura 


afetada. 


21 METODOLOGIA 


A metodologia selecionada para a formulação do trabalho é considerada altamente 
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rica didaticamente. Dessa forma, foram necessárias quatro etapas para a realização do 
trabalho, as quais são apresentadas de forma esquemática na Figura 1. 


Inspeção/Vistoria nas 
estruturas de concreto 


Resultados e 


manifestações Análises p, ) 
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Pesquisa das N e ER cen 
(4 Causas prováveis das 
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A patológicas 


Conclusão e 
Considerações finais 





Figura 1. Fluxograma das etapas executadas no desenvolvimento do trabalho. 


Fonte. Autores (2021) 


Na primeira etapa (Inspeção/Vistoria nas estruturas de concreto) utilizada no 
desenvolvimento deste trabalho, foi iniciada com a visita técnica dos alunos em estruturas 
de concreto armado, neste momento, decidiu-se fazer inspeções/vistoria mais detalhadas 
em estruturas dos elevados (Mariguela; Daniel Berg e Gunnar Vingren) na Região 
Metropolitana de Belém (RMB). Esta etapa consistiu na verificação das manifestações 
patológicas utilizando os sentidos da visão, da audição e do tato. O objetivo foi identificar 
e registrar através de fotografias as manifestações patológicas existentes no sistema 
estrutural de concreto armado. 

Na segunda etapa (Pesquisa das manifestações patologicas) consistiu em realizar 
pesquisas relacionadas a descrição dos aspectos gerais das manifestações patológicas 
encontradas na literatura e apresentadas nos sistemas estruturais em concreto armado; 
acompanhadas da terceira etapa (Levantamento das causas prováveis das manifestações 
patológicas) que ocorreram durante a realização de inspeção/vistoria nas estruturas de 
concreto. 

Na quarta etapa (Resultados e Análises) foram apresentadas as hipóteses de 
causas prováveis das anomalias identificadas nas etapas anteriores com formulações de 
hipóteses baseadas na semelhança dos casos encontrados na literatura e outras citadas 
pelos autores referenciados neste trabalho. Por fim, a quinta etapa acompanha a realização 
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da análise dos resultados e apresenta conclusões do trabalho. 

Ressaltamos que algumas limitações foram encontradas no desenvolvimento 
do trabalho, merecendo destaque as informações incompletas referentes as etapas de 
execução da construção, a impossibilidade de obtenção de amostras através de processo 


destrutivo e a acesso/inexistência de projetos complementares. 


31 RESULTADOS 


O critério de escolha dos elevados objeto de estudo partiu do quesito estrutural, ou 
seja, a definição do concreto armado como material constituinte do sistema estrutural e 
ter suas estruturas aparentes. Pois nestas condições seria possível identificar e analisar, 
a partir das visitas a campo e fotografias digitais, as eventuais manifestações patológicas. 


3.1 Elevados da Região Metropolitana de Belém 


O projeto intitulado “Ação Metrópole” do Governo do Estado do Pará que objetiva 
impantar o Sistema Integrado de trasporte da Região Metropolitana de Belém tem como 
finalidade solucionar o grande problema de mobilidade da Cidade de Belém, uma das 
primeiras etapas do projeto já foi concluida que corresponde a construção do elevado 
“Gunnar Vingren” nas Avenidas Júlio Cezar e Centenário, e “Daniel Berg”, nas Avenidas 
Júlio Cezar e Pedro Álvares Cabral, Figura 2. Além disso, há outros elevados que estão 
no mesmo contexto de mobilidade do sistema integrado de transporte como o elevado 
“Mariguela” localizado nas Avenidas Dr. Freitas e Almirante Barroso; Elevados do “Complexo 
do Entroncamento”, nas Avenidas Almirante Barroso e Rodovias Augusto Montenegro e Br- 
316; e Elevado “Mário Covas”, nas Rodovias BR 316 e Mário Covas, conforme mostra a 


Figura 2. 
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Figura 2. Localização dos Elevados na Região Metropolitana de Belém, da esquerda para a direita. 


Fonte. Autores (2021) 


3.2 Elevado “Gunnar Vingren” 


O elevado “Gunnar Vingren” está localizado nas Avenidas Júlio Cezar e Centenário, 
no bairro Val-de-Cans, Figura 3, se refere a uma ação governamental conhecida como 
Ação Metrópole a fim de compor um sistema integrado de transporte viário na RMB. 


EE — 





Figura 3. Elevado “Gunnar Vingren”. 


Fonte. DOL — Diário on line. (2016) 
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O Elevado em questão foi constituído em estrutura em concreto armado (pilares, 
vigas, lajes maciças e proteção). O local já foi palco de acidentes, inclusive, com registros 
de mortes, onde um condutor de motocicleta caiu do elevado de uma altura de cerca de 6 
metros, após colidir com a mureta de proteção. Segundo o jornal DOL — Diário on line (2016), 
os condutores afirmam que o elevado tem a pista estreita, com sinalização precária, além 
de ter mureta baixa e pouca iluminação pública, o que prejudica ainda mais o tráfego no 
período noturno. A predisposição da estrutura, ou de uma de suas partes, para apresentar 
manifestações patológicas pode ser originada durante a fase de projeto, de construção ou 
ser adquirida na fase de uso. 

Com base nos dados de campo, as manifestações patológicas mais evidentes 
foram provenientes de fissuras, nichos, irregularidades geométricas, manchas escuras, 
exposição e corrosão da armadura, eflorescências e lixiviação. A seguir, no Quadro 1, 
serão apresentadas algumas das anomalias observadas. 


Manifestações Patológicas Imagens obtidas 1n foco 


Fissuras (1) 


Manchas escuras (3) 
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Quadro 1. Registro das manifestações patológicas identificadas no Elevado “Gunnar Vingren”. 


3.3 Elevado “Daniel Berg” 


O elevado “Daniel Berg” está localizado nas Avenidas Júlio Cezar e Pedro Álvares 
Cabral, no bairro Val-de-Cans, Figura 4, se refere a uma ação governamental conhecida 
como Ação Metrópole a fim de compor um sistema integrado de transporte viário na RMB. 


e code 
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Figura 4. Elevado “Daniel Berg”. 


Fonte. DOL — Diário on line (2016). 
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O Elevado em questão foi constituído em estrutura em concreto armado (pilares, 
vigas, lajes maciças e proteção). Com base nos dados de campo, as anomalias mais 
evidentes foram provenientes de choque mecânico, agentes biológicos, fissuras, manchas 
escuras, exposição e corrosão da armadura, eflorescências e lixiviação, conforme mostra 
a Figura 5. 





(A) (B) 
Figura 5. Anomalias. (A) Ação mecânica; (B) Infiltração e corrosão; (C) Manchas Escuras. 


Fonte. Autores (2021) 


3.4 Elevado “Mariguela” 


O elevado “Mariguela” está localizado nas Avenidas Almirante Barroso e Dr. Freitas, 
no bairro Marco. O Elevado em questão foi constituído em estrutura em concreto armado 
(pilares, vigas, lajes maciças e proteção). Com base nos dados de campo, as manifestações 
patológicas mais evidentes foram provenientes de causas naturais, fissuras com infiltrações, 
manchas escuras e exposição e corrosão da armadura, conforme mostra as Figuras 6. 





(A) (B) (C) (D) 


Figura 6. Elevado Mariguela. (A) Fungos; (B) Manchas; (C) Exposição e corrosão de armaduras e (B) 
Ações Mecânicas. 


Fonte. Autores (2021) 
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3.5 Complexo de Elevados do Entroncamento 


O complexo de elevados do entroncamento está localizado nas Avenidas Almirante 
Barroso, Rodovia Augusto Montenegro e Rodovia BR-316, no bairro Souza. Os Elevados 
em questão foram constituídos em estrutura em concreto armado (pilares, vigas, lajes 
maciças e proteção). Com base nos dados de campo, as manifestações patológicas mais 
evidentes foram provenientes de infiltrações, manchas escuras, eflorescências e lixiviação 
e causas naturais, conforme mostra as Figuras 7. 





Figura 8. Complexo de Elevados do Entroncamento. (A) Infiltração; (B) Manchas Escuras; (C) 
Eflorescência e Lixiviação; (D) Fungos. 


Fonte. Autores (2021) 


3.6 Elevado “Mario Covas” 


O elevado “Mario Covas” está localizado na Rodovia BR-316, no município de 
Ananindeua, Km 4, Figura 8. Também conhecido como “viaduto do coqueiro”, possui um 
trevo incompleto com 3 pétalas, inaugurado no início da década de 1980 e, atualmente, 
está em obras com a construção da quarta pétala que interligará a Avenida João Paulo Il. 
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Figura 8. Elevado “Mário Covas”, também conhecido como “Viaduto do Coqueiro”. 


Fonte. DOL — Diário on line. (2016) 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 4 


O Elevado em questão foi constituído em estrutura em concreto armado (pilares, 
vigas, lajes maciças e proteção). Com a investigação, as anomalias mais evidentes estão 
na Figura 9, podendo ser originadas durante a fase de projeto, de construção ou adquiridas 
na fase de uso. 





“—B 
Figura 9. Elevado Mário Covas. (A) Infiltração; (B) Manchas escuras; (C) Fissuras e deterioração. 


Fonte. Autores (2021) 


41 DISCUSSÕES 


Através das inspeções realizadas nos elevados, foi possível observar que as 
manifestações patológicas apresentam fatores variados e diversificadas anomalias como 
corrosão, fissuras, infiltrações, presença de umidade contínua, manchas, fungos, ações 
mecânicas (choque), nichos, lixiviação e falta ou deficiência de manutenção. 


1) Os elementos dos sistemas estruturais estudados (neste caso, pilares, 
vigas e lajes) estão expostos à ação de intempéries, e com isso, a deterioração dos 
materiais que compõem o sistema estrutural reagem com os agentes externos e 
acarretam, muitas vezes, problemas de desgastes variados. Portanto, é concebível 
considerar que estes fatores relacionados a um mau planejamento e/ou execução, 
acabam por contribuir com as principais manifestações patológicas apresentadas 
neste trabalho. 


2) As fissuras e infiltrações foram detectados nos elementos estruturais 
analisados (pilares, vigas e lajes), e são pertinentes pois chamam mais atenção dos 
leigos, condutores e usuários, sendo indicativo de anormalidade no sistema ou que 
algo está para acontecer. A caracterização da fissuração como deficiência estrutural 
dependerá da origem, intensidade e magnitude do quadro de fissuração existente, 
embora o concreto, por ser material com baixa resistência à tração, fissurará por 
natureza, sempre que as tensões trativas, podendo ser instaladas pelos mais 
diversos motivos, superando a sua resistência última à tração (SOUZA e RIPPER, 
1998; MARTINS e FIORITTI, 2016). 
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3) Os nichos de concretagem e segregações no concreto dos elementos 
estruturais estão particularmente ligados à execução e supervisão inadequada 
da concretagem das peças, e apresentam como causas prováveis: a baixa 
trabalhabilidade do concreto, insuficiência no transporte, no lançamento e no 
adensamento do concreto, além de alta densidade de armaduras, corroborando 
com analise de Martins e Fioritti (2016). 


4) Outra anomalia observada foi a desagregação do concreto tanto nos pilares 
quanto nas vigas e lajes, principalmente no Elevado “Mariguela”, esta patologia 
consiste na perda de massa de concreto devido a um ataque químico expansivo 
de produtos inerentes ao concreto e/ou devido à baixa resistência do mesmo, 
caracterizando-se por agregados soltos ou de fácil remoção, causando separação 
de partes do concreto, provocada, em geral, pela expansão devido à oxidação ou 
dilatação das armaduras, e também pelo aumento de volume do concreto quando 
este absorve água (SILVA, 2011: VITÓRIO, 2005 e 2003). 


5) As Manchas escuras foram percebidas nos elementos estruturais e podem 
ser causadas pelos seguintes fatores: deficiência dos dispositivos de drenagem, 
levando a ocorrência de infiltração e umidade, estimulando para que microrganismos 
se proliferem nos mesmos, implicando no surgimento de bolor, sendo agravado 
principalmente pela falta de manutenção nas estruturas. As manchas na superfície 
do concreto alteram a sua textura e uniformidade de coloração causando prejuízos 
estéticos, podendo ainda, gerar o desenvolvimento de problemas patológicos mais 
sérios. As manchas brancas na superfície do concreto das vigas e lajes são definidas 
como eflorescências e lixiviação, respectivamente; e essas alteram a uniformidade 
de coloração do concreto e sua textura, causando prejuízos estéticos. (SILVA, 
2011). O processo de lixiviação ocorre na presença de água quando o concreto for 
mal adensado, e/ou estiver, por algum motivo, fissurado ou apresentar juntas mal 
executadas, permitindo a penetração da água. 


6) A exposição e consequente corrosão de armaduras nos elementos 
estruturais (pilares, vigas, lajes e proteção) tem suas causas prováveis associadas 
principalmente à utilização de um concreto de baixa qualidade, utilização de água 
ou agregados não puros, possíveis erros projetuais bem como erros de execução, 
durante a vibração do concreto e na execução do recobrimento da armadura, assim 
como ação mecânica e uma manutenção precária desse concreto como outros 
problemas relacionados com as intempéries, como excesso de umidade e o contato 
com águas residuais. 


5 | CONCLUSÃO 


Dos resultados apresentados, a manifestação patológica que mais se destacou, 
tanto pela quantidade quanto pela intensidade, foram as fissuras. Suas consequências 
estão relacionadas à durabilidade das estruturas de concreto armado, pois são veículos 
de entrada de água e agentes agressivos para o interior da massa, e com isso, acarretam 
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em patologias mais sérias, como por exemplo, a desagregação do concreto e a corrosão 
de armaduras. 

As manchas escuras também são encontradas em vários fatores analisados, 
ressaltamos que mesmo não apresentando grandes riscos para o bom funcionamento do 
sistema como um todo, podem acarretar no desenvolvimento de problemas patológicos 
mais sérios, ou intensificar aqueles já existentes. 

A exposição e o avançado processo de corrosão de armaduras foi a anomalia que 
mais chamou a atenção, tanto pela quantidade de elementos com esse tipo de manifestação 
patológica, como também pela situação crítica em que alguns desses elementos já se 
encontravam. 

Em menor proporção estão o processo de lixiviação em vigas e lajes, os nichos 
e as irregularidades geométricas, fenômenos que ocorrem normalmente, por dissolução 
e o carreamento do hidróxido de cálcio existente no cimento por meio da ação da água, 
ocorrendo eflorescências de carbonato e redução do pH do concreto; quando o concreto 
foi mal adensado, e/ou esteja, por algum motivo, fissurado ou apresente juntas mal 
executadas, permitindo a penetração da água. Estas manifestações patológicas precisam 
da mesma forma como as abordadas anteriormente, de cuidados e recuperação para não 
se agravarem com o passar do tempo. 

Portanto, com o estudo e as análises realizadas, é possível afirmar que muitas 
das manifestações patológicas detectadas poderiam ter sido minimizadas, e até mesmo 
evitadas, caso houvesse maior controle das etapas do processo construtivo, e que ao 
serem associadas a uma efetiva manutenção das estruturas de concreto armado, poderiam 
proporcionar maior qualidade e vida útil aos elementos estruturais. À falta de manutenções 
adequadas as estruturas, que não são realizadas de forma constante e preventiva, 
contribuindo para o aparecimento das anomalias, como também para o aumento do grau e 
intensidade delas, acarretando na necessidade de processos de reparos e reestruturação, 
afetando o funcionamento da edificação e levando a um gasto financeiro maior do que se 


houvesse uma prevenção efetiva e constante. 
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RESUMO: Neste capítulo descreve-se a 
formulação co-rotacional de um elemento de 
viga unificado que engloba as teorias de Euler- 
Bernoulli e de Timoshenko e que não apresenta 
bloqueio por cisalhamento. A cinemática co- 
rotacional se baseia na separação do movimento 
de um sólido em uma parte deformacional, 
e a outra, em movimento de corpo rígido. O 
movimento deformacional do elemento de viga é 
descrito por três modos naturais de deformação 
relacionados aos esforços axial, de flexão pura e 
de flexão simples, respectivamente. Os esforços 
internos gerados pelos modos de deformações 
naturais são auto-equilibrados, o que permite 
obter uma matriz de rigidez tangente consistente. 
Descreve-se de maneira detalhada a obtenção 
das matrizes de rigidez material e geométrica. 
Aplica-se o elemento de viga unificado na 


Coleção desafios das engenharias: Engenharia civil 4 


análise não linear de arcos, demonstrando que 
esse elemento é capaz de lidar com grandes 
translações e rotações de corpo rígido. 
PALAVRAS-CHAVE: Formulação co-rotacional; 
Elemento de viga Bernoulli-Timoshenko; Modos 
de deformação naturais; Métodos dos Elementos 
Finitos; Não-linearidade geométrica. 


NONLINEAR ANALYSIS OF ARCHES 
USING CO-ROTATIONAL FORMULATION 


ABSTRACT: This chapter describes the co- 
rotational formulation of a unified beam element 
that combine Euler-Bernoulli and Timoshenko 
theories, which don't have problems with shear 
locking. The co-rotational kinematics is based 
on separation of the montion in deformational 
and rigid body componenis. The deformation 
movement of the element is described by three 
natural modes of deformation related to the 
axial efforis, pure bending and simple bending. 
The internal forces generated by the natural 
deformation modes are self-balanced, which 
allows to obtain a consistent tangent stifífness 
matrix. Developmentofthe material and geometric 
stiffness matrix is described in details. The unified 
beam element is applied in the nonlinear analysis 
of arches, demonstrating that it has the ability to 
handle large translations and rigid body rotations. 
KEYWORDS: Co-rotational formulation; 
Bernoulli-Timoshenko element beam; Natural 
modes of deformation; Finite Element Methods; 
Geometric nonlinearity. 


Capítulo 5 


11 INTRODUÇÃO 


A ideia principal da formulação co-rotacional é a decomposição do movimento de 
um sólido na soma de um movimento de corpo rígido e de um movimento deformacional. 
Esta formulação é um caso particular da descrição Lagrangiana para levar em conta a 
não linearidade geométrica na análise estrutural fazendo o uso do método dos elementos 
finitos. Assume-se nesta formulação que as translações e rotações de corpo rígido poderão 
ter qualquer magnitude porém as deformações deverão ser pequenas. Por esta razão é 
possível utilizar elementos finitos lineares em problemas que envolvam a não linearidade 
geométrica, sendo esta a principal motivação para o uso da formulação co-rotacional. 
Nos anos 70 e 80 esta formulação foi utilizada na análise não linear de diversos sistemas 
estruturais utilizados na engenharia aeronáutica e aeroespacial, na engenharia naval e 
mecânica, entretanto tendo pouca aplicação na engenharia civil nessa época. Em meados 
dos anos 90 esta formulação foi utilizada na análise não linear física, bem como, na 
análise dinâmica não linear de diversos sistemas estruturais. Skallerud et al. [1] utilizaram 
a formulação co-rotacional para analisar o comportamento de dutos submarinos expostos 
a grandes mudanças de geometria bem como a falhas locais. Yaw et al. [2] aplicaram a 
formulação co-rotacional nos métodos sem malhas para analisar sólidos 2D sujeitos as 
não linearidades geométrica e física. Felippa e Haugen [3] descrevem uma formulação co- 
rotacional unificada para elementos finitos de vigas, placas e cascas. Por fim, podem ser 
citados outros autores [4,5,6,7] que utilizaram a formulação co-rotacional para análises não 
lineares de diversos sistemas estruturais. Neste trabalho apresenta-se a formulação co- 
rotacional de um elemento de viga 2D com as seguintes particularidades: a) Define-se um 
vetor de deslocamentos de corpo rígido através de um sistema de referência fixo no centroide 
do elemento; b) Este elemento pode sofrer rotações de corpo rígido de qualquer magnitude 
devido ao uso da função trigonométrica módulo; c) Adotam-se três modos de deformações 
naturais [8,9] para representar a parte deformacional do movimento do elemento; d) 
Estes modos de deformações naturais representam esforços internos auto-equilibrados; 
e) Obtém-se o vetor de forças nodais a partir dos esforços internos auto-equilibrados; f) 
Utiliza-se o princípio dos trabalhos virtuais complementar [10,11] para obter a matriz de 
rigidez elástica que integra as teorias de vigas de Euler-Bernoulli e de Timoshenko; 9) Este 
elemento será aqui denominado de elemento de viga 2D unificado; h) Este elemento não 
apresenta bloqueio por cisalhamento; i) Descreve-se detalhadamente o desenvolvimento 
algébrico para o cálculo das matrizes de rigidez co-rotacional e geométrica, assim como 
a expressão completa da matriz de rigidez tangente. Finalmente, se aplica este elemento 
para a análise não linear de arcos planos. 
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21 DESCRIÇÃO CO-ROTACIONAL 


Seja um sistema de coordenadas globais cuja base é ortonormal conforme mostra- 
se na Figura 1. 





Ec 


Figura 1: Movimento do elemento de viga 2D. 


Para expressar as variáveis cinemáticas na configuração indeformada são utilizadas 
as coordenadas materiais (X, Y). Para expressar as variáveis cinemáticas na configuração 
deformada são utilizadas as coordenadas espaciais (x, y). Na configuração indeformada 
as coordenadas nodais do elemento são dadas por (X, X,) e (X,, X,), respectivamente. A 
longitude e inclinação inicial são dadas por: 


lo = JX +Yã 


lo — Ao: 
Po = Zarctan 2) 
ha 


(1) 


— emqueX, =X,-X e Y, = Y,-Y,. Para obter a expressão (1b) foi utilizada a relação 
tan (- Po). Esta expressão é singular se Y,=0. Se X,=l, assume-se y,=0. Se X, =-|, 
assume-se py =It De maneira similar, na configuração deformada as coordenadas nodais 
do elemento são dadas por (x, Y,) e (X,, Y,), respectivamente. A longitude e inclinação 
atuais são dadas por: 
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[= JX% + Yã 


[— X94 
p = Zarctan —+ 
Y21 


em que x, =X,-X ey, =Y,-Y,.. Para obter a expressão (2b) foi utilizada a relação 


(2) 


a | | 
tan (- Po). Esta expressão é singular se y, =0. Se x, =|, assume-se qp =0. Se x, =-l, assume- 
se p=r. Como mostra-se na Figura 1, o movimento do elemento de viga 2D é dado pelo 
vetor de deslocamentos nodais, expresso em coordenadas globais, como: 


u" = (u,v,0,u20202) (3) 


As coordenadas nodais na configuração deformada podem ser escritas em 
função dos deslocamentos nodais como x =X +u e x,=X,+u,, respectivamente. A partir 
destas expressões obtém-se que x, =X, +U, € ),=Y,+Y,, em que u,=U,-U € V,=V,vV.,, 
respectivamente, o que permite reescrever que -vV(X,, + u5,)* + (1, + v5,)*. Para 
definir o movimento de corpo rígido é fixado no centroide do elemento de viga 2D um 
sistema de coordenadas locais (x, y.) que o acompanha durante todo o movimento, 
conforme mostra-se na Figura 1. Nota-se que o movimento de corpo rígido é definido por 
duas translações e uma rotação do sistema de referência local entre as configurações 


deformada e indeformada, e que se escreve como: 


1 
2 (u, + ui) | 


ft, = 1 (4) 
2 Uz + Vi) 


P— Po 
A rotação de corpo rígido entre as configurações deformada e indeformada é dada 


por y-q,. Por outro lado, o movimento deformacional em coordenadas locais, conforme 


mostra-se na Figura 1, é definido como: 
Ua — 04 (5) 


em que 9,=0,-(p-p,) e 9,=0,(-,) são rotações deformacionais dos nós 1 e 2, 
respectivamente, e u=/-|, é o estiramento ou encurtamento do elemento. Posteriormente, 
serão adotadas as rotações deformacionais geradas pelos esforços de flexão pura e 
simples, respectivamente. 


2.1 Modos de deformações naturais 


Foram necessários três graus de liberdade para definir o movimento de corpo rígido 
do elemento de viga 2D de acordo com o vetor de deslocamentos nodais dado pela equação 
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(4). Portanto, os três graus de liberdade restante definem o movimento deformacional do 
elemento de acordo com o vetor de deslocamentos nodais dado pela equação (5). Neste 
trabalho, adota-se o conceito dos modos de deformações naturais definidos nas referências 
[8,9] como mostra-se na Figura 2. O primeiro modo é obtido pela aplicação de um esforço 
axial nos nós do elemento conforme mostra-se na Figura 2a. O segundo modo é obtido ao 
considerar o elemento sujeito a um estado de flexão pura o que implica na ausência de 
esforço cortante. Este modo define uma linha elástica simétrica em relação ao eixo local y, 
conforme mostra-se na Figura 2b. No terceiro modo o elemento está sujeito a um estado de 
flexão simples que resulta em um esforço cortante constante. Este modo define uma linha 
elástica antissimétrica em relação ao eixo local conforme mostra-se na Figura 2c. 





O u 
'L 


Figura 2: Modos de deformações naturais. a) estiramento. b) elástica simétrica. 
c) elástica antissimétrica 


Portanto, pode-se reescrever o vetor de deslocamentos deformacionais dado em 


(5), em função dos modos de deformações naturais, da seguinte maneira: 


ui = (u,0,,04) (6) 
em que: 
0, =0,-—04 (7) 


da =0,+0,— 2(P — Po) 


Para obter as expressões (7b) e (7c) utilizam-se as rotações deformacionais 0, e 
9, dadas em (5). Logo, conforme as Figuras 2b e 2c, somam-se a rotações simétricas e 
antissimétricas de cada nó do elemento. Por último, igualam-se essas somas às rotações 
0, e 6,, respectivamente. Este procedimento resulta no seguinte sistema de equações: 
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] di 
—58s + 5% E 0, E (gp — Po) 
; (8) 
2 Os +50 =0—(P—Yo) 


cuja solução obtém os valores de 0, e 6 dados pelas equações (7b) e (7c), 
respectivamente. Para que não haja restrições quanto as magnitudes dos ângulos, a 
expressão da rotação antissimétrica deve ser calculada usando a função trigonométrica 


módulo [11], tal que: 
da = laços (7 e T) =X (9) 


A função módulo coloca um valor arbitrário no intervalo [0,27], e o último termo 
restabelece a simetria em relação ao valor zero. A omissão deste passo pode levar a 
problemas de divergência quando o elemento de viga 2D sofre rotações, trT, *271, *3TT.... 

Na Figura 3 mostra-se a aplicação da função módulo dada em (9) para ângulos que 
variam de zero até 6=20r1, que corresponde a 10 voltas em círculo de raio unitário. Nota-se 
que os valores estão dentro do intervalo -r1 <0*<rr devido à simetria em relação ao valor 


rr 


Figura 3: Função módulo para ângulos arbitrários. 


Zero. 


9º=mod(0 +7) x 


2.2 Vetor de forças internas 


É importante ressaltar que os modos de deformações naturais geram esforços auto- 


equilibrados como se mostra na Figura 4. O vetor de esforços auto-equilibrados é expresso 
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como: 


fa = (nmma) (10) 


Ne 


» a 


Figura 4 — Esforços auto-equilibrados. a) Esforço axial. b) Flexão pura. c) Flexão simples. 


fu 








Figura 5 — Vetor de forças nodais. a) Em coordenadas globais. b) Em coordenadas locais. c) Auto- 
equilibradas. 


A combinação dos esforços auto-equilibrados, mostrados na Figura 4, gera o vetor 
de forças nodais em coordenadas locais como mostram as Figura 5b e 5c. O vetor de forças 
nodais também pode ser expresso em coordenadas globais como mostra a Figura 5a. A 
relação entre o vetor de esforços auto-equilibrados e o vetor de forças nodais expresso em 


coordenadas locais é dada por: 


fé = Sfa (11) 
em que S é a matriz de transformação dada pela seguinte expressão: 
—1 0 O 
Oo 0 2/4 
midptto Sd d 4 
S = 10 1 (12) 
O 0 -2/ 
0 1 1 


Analogamente, a relação entre o vetor de forças nodais em coordenadas locais, 
ver Figura 4b, e o vetor de forças nodais em coordenadas globais, ver Figura 4a, é escrita 


como: 
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f=Rf (13) 
em que R é a matriz de rotação do sistema de coordenadas locais para o sistema de 


coordenadas globais. Define-se esta matriz através da seguinte expressão: 


cosp -senp O 0 0 0 

senp cosp O OQ 0 0 
R=| 0 0 1 0 0 0 (14) 

0) 0 0 cosy -senp 0 

O 0 0 senp cosy 0 

O 0 0 0 0 1 


Finalmente, substituindo a equação (11) na equação (13), obtém-se a relação entre 
o vetor de forças auto-equilibradas e o vetor de forças nodais em coordenadas globais, tal 
que: 
f = R$fa (15) 
Por outro lado, os deslocamentos nodais podem ser escritos em relação ao sistema 


de coordenadas locais como: 


ut =R'u (16) 


31 MATRIZ DE RIGIDEZ TANGENTE 


A relação entre os incrementos de deslocamentos e de forças em coordenadas 


locais se expressa como: 
df* = Kºduº (17) 
em que Kº é a matriz de rigidez tangente. Utilizando as equações (13), (16) e (17), 
determina-se a matriz de rigidez tangente em coordenadas globais da seguinte maneira 


K = RKºR' (18) 
Ao aplicar um incremento infinitesimal de deslocamentos a partir da configuração de 
equilíbrio atual conforme mostra-se na Figura 6, os incrementos de rotação de corpo rígido 
e de longitude do elemento se escrevem como: 


— dvj-dvy 
Sp 


dy dl = dus — duí (19) 


em que (duj, dvj, dus, dv;) são as componentes dos deslocamentos nodais em 


relação ao sistema de coordenadas locais. 
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Figura 6 — Movimento incremental do elemento viga 2D. 


Considerando que o trabalho interno independe do sistema de coordenadas, e 


levando em conta a equação (11), a variação do trabalho interno pode ser escrito como: 
Fr Tr 
dV = du” fº = du Sfa = duify (20) 


em que dui=(du, do, do) é a variação dos modos de deformações naturais. A 
última igualdade na equação acima se cumpre para qualquer valor arbitrário de f,, como 


consequência se obtém que: 


A variação do vetor de forças internas dado pela equação (15) se escreve como: 
df = RSdf, + RdSfy + dRSf, (22) 
com, 
-senp -—cosp 0 0) 0 0 O Q O | 
cosp -senp O 0 0 0 0 O -2/1º 
E O 0 0 0 0 O ao 0 0 29 
did O Õ O —senp -—cosp O sai dd 0 0 [Ê tea) 
0 0 O cosp -senp 0 O 0 2H 
O O O 0 O O 00 0 


A relação entre as variações dos esforços internos e dos modos de deformações 


naturais se expressa como: 
dfa = Kadu, (24) 


em que K, é a matriz de rigidez elástica, cuja dimensão é 3 x 3. Levando em conta 
as equações (13), (21) e (24), a equação (22) pode ser reescrita como: 


dfº = SKuS” du" + (dS + RT dRS)f, (25) 


Finalmente, levando em conta a equação (17), a matriz de rigidez tangente em 
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relação ao sistema de coordenadas locais pode ser escrita como: 
Ke =SKiST +K, (26) 


O primeiro termo leva em conta os coeficientes de rigidez devido aos modos de 
deformações naturais, enquanto que K é a matriz de rigidez co-rotacional, que representa 
o efeito da variação da rotação de corpo rígido do sistema de referência local y do efeito do 
esforço cortante devido à variação da longitude do elemento de viga 2D. De acordo com a 
equação (25), a matriz de rigidez co-rotacional é definida como K =dS+R'dRS. De acordo 
com a equação (23) verifica-se que a variação da matriz de transformação S depende da 
variável dl e que a variação da matriz de rotação R depende de dy. Levando em conta a 


equação (19), depois de um desenvolvimento algébrico, chega-se a: 


O Q O 0 -Q M 

7 Q N O-Q -N 0 
= a) O O 0 O 0 27 
*=Tlo Q00 00 dá 

=) -N O Q N 0 


- O O 00 O Ol 
Esta matriz forma parte da matriz de rigidez geométrica completa do elemento de 


viga 2D como se demonstrará na seção 3.8. 


3.1 Matriz de rigidez elástica 


O primeiro modo de deformação é o alongamento/encurtamento do elemento de 
viga 2D e seu coeficiente de rigidez, de acordo com a teoria da resistência dos materiais, 
se obtém através da relação incremental dN= du, Para o elemento sujeito a flexão pura 
obtém-se o coeficiente de rigidez através do Princípio dos Trabalhos Virtuais Complementar 
(PTVC), o que resulta na seguinte relação incremental dM, = 6, O terceiro modo de 
deformação natural representa o elemento sujeito a flexão simples o que resulta em um 
esforço cortante constante. Desta maneira, ao aplicar o PTVC, obtém-se que: 


l 


MG? à 1M2l QU 
Mat = | ( pas a q O (28) 
aZa EL "GA, > EL CA 





O 
em que A, é a área da seção transversal corrigida pelo fator de forma que leva em 
conta o efeito da distribuição das tensões de cisalhamento na seção transversal. Introduz- 


se o efeito da deformação por cortante através dos seguintes parâmetros: 


x 1 ias 12E!I (29) 
da = Trop — GAo!2 
Usando a expressão do esforço cortante qa e a equação (28), obtém-se a 


seguinte relação incremental: 
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El 
dMa = 3a 7 dO (30) 


Finalmente, a matriz elástica é definida como: 


EA 
— 0 0 
EI 
K=|0 — 0 (31) 
El 
O O 3a 


Cabe destacar que os coeficientes de rigidez a flexão desta matriz levam em conta 
tanto a hipótese de Euler-Bernoulli quanto a hipótese de Timoshenko, por esta razão a 
denominação de elemento de viga unificado. 


3.2 Matriz de rigidez geométrica 


Ademais da matriz de rigidez elástica é necessário levar em consideração a matriz 
de rigidez geométrica devido ao efeito do esforço axial quando o elemento de viga 2D sofre 
um incremento de rotação de corpo rígido. Neste trabalho, para obter os coeficientes da 
matriz de rigidez geométrica considera-se a equação diferencial da viga-coluna sem cargas 
transversais e não levando em conta o efeito do esforço cortante. Considera-se o elemento 
de viga simplesmente apoiado para que se cumpra as condições de contorno do segundo 
e terceiro modos de deformações naturais como se mostra nas Figuras 2b e 2c. Desta 


maneira, a equação diferencial que governa o problema se escreve como: 
ElvC — Nv" =0 (32) 
Após aplicar duas vezes a integração por partes, obtém-se a seguinte forma integral 
desta equação diferencial: 
l 
[over + óv'Nvdx = [óvQ + 69M], (33) 
0 
A relação incremental desta forma integral, após o desenvolvimento algébrico, se 
expressa como: 
| 
| cowerar” + óv Ndv')dx = [óvdQ + ô0dM!]s, (34) 
0 


Para calcular o coeficiente de rigidez geométrico devido ao segundo modo de 


deformação natural assume-se a seguinte igualdade: 
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| 
) ôv'Ndv'dx = [89,dMJ! (35) 
Ú 


Interpola-se o campo de deslocamento incremental dv, o campo de deslocamento 
virtual ôv e suas derivadas de primeira ordem utilizando a função de forma para a elástica 


simétrica definida na Figura 7. 


Elástica simétrica Elástica antissimétrica 
l 
| Vo TEN 
v= -a( baia EO, v 9 (1 é )é8, 
Regra da cadeia 
dC) dC) dg 


dx — de dx 


Figura 7: Função de forma. 


1 1 
Em seguida aplica-se a regra da cadeia obtendo que av=" 606, e ôv= 7500, 


Posteriormente, integrando a equação (35) no espaço da coordenada é chega-se a: 


Il l | 
dM, — 2 Nidda — ka — aNi (36) 


Para calcular o coeficiente de rigidez geométrico devido ao terceiro modo de 


deformação natural, impõe-se a seguinte igualdade: 
| 
| ôv'Ndv'dx = [66,dM,|]! (37) 
O 


Seguindo os passos anteriores para o caso da FSito antissimétrica, ao aplicar 
a regra da cadeia obtém-se que Er IRS] do e id da ALE Desta maneira, 
integrando a equação (37) no espaço da coordenada &€ chega-se a seguinte expressão: 
dM > Ntdo k NI (38) 
= — => — — — tod 
22 20 é o | 


Portanto, define-se a matriz de rigidez geométrica como: 
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(39) 


3.3 Matriz de rigidez tangente completa 

Para obter a matriz de rigidez tangente completa do elemento de viga 2D estendem- 
se as dimensões das matrizes K, e K, para suas formas completas, isto é, para matrizes 
de ordem 6x6 utilizando-se a matriz de transformação S. Consequentemente, a forma 


completa da matriz de rigidez tangente se escreve como: 


Kº =SKgST + SK9ST + K,=K$+K$ (40) 
com, 
EA o : EA 
; 
EI EI EI 
0 12h E Paz 0 1a 6a 77 
El El EI EI 
ds ' Sa [£ Cia + 7 0 ova [2 (Va a DE 
Rd=] EA EA 
0 12Wa 7ã bar O ya az 
f| EI El EI 
O Cha Gha-D7 0 Sha Chat DT] (41) 
0 Q 0 0 -Q 0 
NO NL on Lui 
e 5 10 Q 5 10 
O E NI 2 NI? 0 ! NI! À NI? 
KE=SKyST + K, => 10 > 10 30 
s=Hd *ÉrZio —Q 0 0 Q 0 
Q EN DONE Q ÉN NI 
| 5 10 5 10 
10 30 10 15 


em que Kº é a matriz de rigidez tangente do elemento de viga 2D em coordenadas 
locais, Ki é a matriz de rigidez material e K4 é a matriz de rigidez geométrica. Cabe destacar 
que Turner et al. [12] foram os primeiros autores a obter a matriz de rigidez material para 
um elemento de viga de Timoshenko sem bloqueio por cortante e que em Przemieniecki 
[13] descreve-se a matriz de rigidez material em função do coeficiente ya, definido em 
(29), através da solução de uma equação diferencial ordinária levando em consideração as 
hipóteses de Timoshenko. Por outro lado, assumindo Q=0 na matriz de rigidez geométrica 
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em (41b) chega-se a matriz de rigidez geométrica padrão do elemento de viga de Euler- 


Bernoulli que se escreve como: 


0 0 O «0.0 0 
O 36 3 0 -36 3l 
Wi SP 42d =. =P 
É — — 
Kg 301] O 0 O 0 0 0 (42) 
-0 -36 -31 0 36 -3 
O 3 -E WO -g 4 


Esta matriz aparece pela primeira vez no artigo publicado por Martin [14] nos anos 60. 


4 | EXEMPLOS NUMÉRICOS 


Para realizar a análise não linear geométrica dos exemplos apresentados nesta 
seção, foi utilizado um programa escrito em linguagem Fortran90 denominado co. 


rotating 2Dbeam.f90 escrito pelo primeiro autor deste artigo. 


Deslocamentos em coordenadas globais 


u” = (u,, 04,04, Us, Vo, 05) 


Configuração deformada 


T=X+0=(X +u,Y +0,0,4 + Xv, +Y,,0,) 


Rotação de corpo rígido: p =p — py, 

P, = 2 arctan (ça : y = 2arctan (Em ) (1), (2) 

21 91 

Modos de deformações naturais 

u=l— lo 

0,=0,-0, 

0,=0,+0, — 2(/ — Po) 
Função módulo para ângulos arbitrários 

= mod(9, +mT)—7 (9) 
Esforços internos 

f,=K,u, (10), (31) 
Esforços em coordenadas globais 

f= RST, (12), (14), (15) 
Matriz de rigidez tangente em coordenadas globais 


K=RKR (14), (18), (40) 





Tabela 1: Algorítmo da descrição co-rotacional. 
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A Tabela 1 resume os principais passos a serem implementados em um programa 
computacional para o cálculo do vetor de forças internas e da matriz de rigidez tangente do 
elemento de viga 2D, ambos expressos em coordenadas globais. 

Em todos os exemplos apresentados nesta seção foi adotada uma tolerância para 
convergência de 10%. O objetivo nesta seção é mostrar o desempenho do elemento de 
viga unificado desenvolvido neste trabalho para reproduzir o comportamento fortemente 
não linear de alguns sistemas estruturais. É importante destacar, como será visto nos 
exemplos a seguir, que o elemento de viga unificado aqui descrito não apresenta bloqueio 
por cortante, além da habilidade do mesmo em lidar com grandes rotações de corpo rígido. 


4.1 Arco circular de grande altura rotulado e engastado 


Como mostra a Figura 8a, este exemplo consiste em um arco de grande altura de 
diretriz circular com uma extremidade engastada e a outra apoiada com um vínculo de 
segundo gênero. Este arco está sujeito a uma carga concentrada no vértice. Na Figura 8a 
detalham-se as propriedades geométricas e mecânicas do arco que foi discretizado com 20 
elementos de viga 2D unificado. 


a) Pa b) 







E = INKH tt : 7 
EB WApNEI 


Ft = lh) 
[0] 215" 





250 UU 


Figura 8: Arco circular de grande altura. a) Propriedades geométricas e mecânicas. b) Trajetórias de 
equilíbrio. 


Na Figura 8b mostram-se as trajetórias de equilíbrio que descrevem a evolução 
dos deslocamentos horizontal u e vertical v do vértice do arco. Foram comparados os 
resultados da trajetória de equilíbrio Px v com os resultados obtidos por Wagner [15], onde 
pode-se observar uma boa concordância entre ambos resultados. Para obter as trajetórias 
de equilíbrio foi utilizado o método de longitude de arco com restrições das iterações em 
trajetória ortogonal nas fases corretoras. A longitude de arco adotada foi de 12 para 100 
passos de carga e 436 iterações totais. O número médio de iterações por passo de carga 
foi de 4,36. Pode-se observar que as trajetórias de equilíbrio mostradas na Figura 8b 
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apresentam pontos limites e turnning points. 


4.2 Arco semicircular biarticulado 


Neste exemplo considera-se um arco semicircular biapoiado sujeito a uma carga P 
excêntrica em relação ao seu vértice como mostra a Figura 9a. Nesta figura encontram-se 
detalhadas as propriedades geométricas e mecânicas do arco que foi discretizado com 50 
elementos de viga unificado (EVU). Na Figura 9b mostra-se a curva da carga P versus o 
deslocamento v do vértice. Observa-se nesta figura que há 9 pontos limites. 


P(lb) 





E= 2000 psi, G= 1000 psi 
I=lin', A= Win” — 40) 
E =1IWin. f=5/6 


— [5] | 





Figura 9: Arco semicircular biarticulado. a) Propriedades geométricas e mecânicas. b) Trajetórias de 
equilíbrio. 


Os valores da carga P nos nove pontos limites foram comparados com os valores 
apresentados por Yang e Kuo [16]. Estes autores utilizaram a formulação Lagrageana 
Atualizada, um elemento de viga que inclui termos não lineares da deformação axial de 
Green-Lagrange e uma malha com 26 elementos de viga 2D. Os valores destas cargas 
estão detalhados na Tabela 2. Nota-se que a maior diferença foi de 3,82% observada no 
ponto limite 9 (LP9). Para obter a trajetória de equilíbrio foi utilizado o método de controle 
variável dos deslocamentos. O valor inicial da longitude de arco foi de 2,2. Foram utilizados 
859 passos de carga com 91 cortes automáticos devido a divergência no processo iterativo 
da fase corretora de alguns passos de carga. O número médio de iterações foi de 4,28. 
Pode-se notar que a trajetória de equilíbrio mostrada na Figura 9b é fortemente não linear 
com a presença de vários limit points, turning points e loopings. 
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PL1 PL2 PL3 PL4 PLS PL6 PL7 PL8 PL9 


Yang y Kuo 5.813 -8.498 16.149 -22.162 38.566 -49.896 64.875 -82.420 104.611 
EVU 9.811 -8.495 16.204 -22.086 38.932 -50.206 66.786 -83.138 108.61 
diferença (%) 0.03 0.04 0.34 0.34 0.95 0.62 2.99 0.87 3.82 


Tabela 2: Valores da carga P nos pontos limites — (lb) 


4.3 Arco circular biarticulado abatido sob carga excêntrica 


Neste exemplo considera-se um arco abatido de diretriz circular com ambas 
extremidades apoiadas e sujeito a uma carga concentrada P com uma excentricidade de 
0,2m em relação ao seu vértice. Desta maneira, considerou-se a carga P e um momento 
fletor M=0,2P aplicados no vértice do arco. Na Figura 10a detalham-se as propriedades 


geométricas e mecânicas do arco que foi discretizado com 20 elementos de viga unificado. 


| 


E = 5.5556:10 N /m” 
AN Pe G = 2.1368.10 N /m” 
— [= 0.0036mº 
A = 0.36mº 
b = 0.6V/3m 
Fy = 0.24/3m 
f=5/6 


e = U.Zm 





P (kN) P (kN) 


LO) bh) | 


o 
7 *ochan e chi 
D / 
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——e 
- « à= 
11 02 04 06 OS 10 12 14 16 0.012-0.008-0.004771 0.004 0.008 0.012 
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Figura 10: Arco circular biarticulado abatido. a) Propriedades geométricas e mecânicas. b) 
Deslocamento vertical do vértice. c) Deslocamento horizontal do vértice. 


Na Figura 10b mostra-se a curva carga P versus deslocamento vertical v do vértice. 
Na Figura 10c mostra-se a curva carga P versus o deslocamento horizontal u do vértice. 
Foram comparados os resultados de ambas trajetórias de equilíbrio com os resultados 


obtidos por Chan e Chui [17] onde pode-se observar uma boa concordância, com exceção 
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de pequenas discrepâncias no tamanho dos laços na direção do eixo x apresentados na 
Figura 10c. Essas discrepâncias, provavelmente, devem-se ao processo de obtenção 
desses pontos. Esses pontos foram obtidos após a digitalização dos gráficos apresentados 
em [17], e também, há de se ressaltar o ajuste das escalas dos eixos coordenados. Para 
obter as trajetórias de equilíbrio foi utilizado o método de controle variável dos deslocamentos 
apresentados na referência [18]. O valor da longitude de arco foi de 0,0085 para 315 passos 
de carga. O número médio de iterações foi de 3,6. Nota-se que as trajetórias de equilíbrio 
do arco mostradas nas Figuras 10b e 10c são fortemente não lineares com a presença de 


alguns limit poinis, turning points e loopings. 


4.4 Arco circular biarticulado abatido sob carga distribuida 


Neste exemplo considera-se um arco abatido de diretriz circular biapoiado e 
submetido a uma carga q uniformemente distribuída ao longo da metade de seu vão, 
conforme mostra-se na Figura 11a. 
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h=2.5m, ] LO mi 
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Figura 11: Arco circular biarticulado abatido. a) Propriedades geométricas e mecânicas. b) Trajetórias 
de equilíbrio. 


Nesta Figura 11a detalham-se ainda as propriedades geométricas e mecânicas do 
arco que foi discretizado com 20 elementos de viga unificado (EVU). Na Figura 11b mostra- 
se a curva da carga normalizada gr/El versus o deslocamento vertical normalizado v/r 
do vértice. Pode-se observar que há quatro pontos limites. Os valores normalizados das 
cargas nos pontos limites foram comparados com os valores obtidos por Xu e Mirmiran [19]. 
Estes autores utilizaram a formulação co-rotacional e um elemento de viga que inclui os 
termos não lineares da deformação axial de Green-Lagrange. Além disso, utilizaram uma 
malha com 20 elementos. Os valores destas cargas estão detalhadas na Tabela 3, onde 
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observou-se que a maior diferença foi de 3,1% no ponto limite 2 (PL2). 


PLI PLZ PLS 
Xu y Mirmiran 13.77 -20.09 33.99 
EVU 13.92 -20.71 3486 





* diferença (%) 109 309 2.56 


4 = " 
Tabela 2: Valores da carga — nos pontos limites 


Para obter a trajetória de equilíbrio foi utilizado o método de controle variável dos 
deslocamentos. O valor inicial da longitude de arco foi de 0,1. Foram utilizados 177 passos 
de carga com 12 cortes automáticos devido a divergência do processo iterativo na fase 
corretora de alguns passos de carga. O número médio de iterações foi de 4,18. Pode-se 
notar que a trajetória de equilíbrio do arco mostrado na Figura 11b apresenta fortes não 
linearidades com a presença de alguns limit points, turning points e um looping. 


5 | CONCLUSÕES 


Neste trabalho mostrou-se, através dos exemplos numéricos, a habilidade do 
elemento de viga unificado em lidar com grandes translações e rotações de corpo rígido. 
Desta maneira, não se produz um aumento artificial da rigidez no comportamento não linear 
de diferentes tipologias estruturais. Foi descrita de maneira detalhada a relação entre a 
formulação co-rotacional e os modos de deformações naturais. A grande vantagem do uso 
dos modos deformacionais naturais é que os mesmos definem a priori esforços internos 
auto-equilibrados, o que permite obter o vetor de forças nodais e a matriz de rigidez 
tangente de maneira consistente, descritos em coordenadas globais, através de operações 
algébricas bastante simples. Ademais, neste trabalho descreveu-se de maneira explícita o 
movimento de corpo rígido do elemento de viga unificado e a dedução da matriz de rigidez 
co-rotacional. Cabe destacar que a principal vantagem da formulação co-rotacional é o 
desacoplamento entre os efeitos locais e globais, e isto permite utilizar uma biblioteca de 
elementos finitos oriundos da análise linear. A utilização de elementos finitos lineares na 
análise não linear geométrica através da formulação co-rotacional se deve ao fato desta 
formulação descrever os efeitos globais como movimento de corpo rígido. Por último, é 
importante ressaltar que distintas fontes de não linearidades físicas podem ser tratadas 


como efeitos locais na formulação co-rotacional. 
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ABSTRACT: The negative impact of inadequate 
final disposal of waste on the environment and on 
health is becoming more and more evident in our 
country. The population has become aware of the 
importance of this problem and has demanded 
more decisive action from public institutions 
to solve it. The alternatives that are being 
offered to solve this problem start today with a 
comprehensive management approach and put 
a lot of emphasis on the sustainability of the 
solutions within legal instruments. Final disposal 
of waste adopts particular characteristics in small 
towns and in rural areas, due to factors such as: 
the lack of resources, by the almost generalized 
subsidy of the cleaning service; the absence of 
information on the negative consequences of 
landfills; ignorance of joint solutions which reduce 
implementation and operation costs of manual 
fillings; the appropriate technology to dispose 
the waste without incurring into greater costs of 
investment and operation; and the absence of 
knowledge about how to deal with the problem of 
inadequate final disposal of waste. 
KEYWORDS: Soil, design, 
management plan. 


environmental 
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VALENZUELA 


11 INTRODUCCIÓN 


Es un hecho cada vez más claro 
en nuestro país el impacto negativo de la 
inadecuada disposición final de residuos en el 
ambiente y en la salud. La población ha tomado 
conciencia de la importancia de este problema y 
ha demandado una acción más decidida de las 
instituciones públicas para solucionarlo. 

Las alternativas que se vienen 


ofreciendo para resolver este problema 
parten hoy de un enfoque de gestión integral 
y ponen mucho énfasis en la sostenibilidad 
de las soluciones y enmarcadas dentro de los 
instrumentos legales. la disposición final de 
residuos adopta características particulares 
en localidades pequehas y en zonas rurales, 
debido a factores como: la falta de recursos, 
por el subsidio casi generalizado del servicio de 
limpieza; la ausencia de información sobre las 
consecuencias negativas de los vertederos; el 
desconocimiento de soluciones conjuntas, que 
reducen costos de implementación y operación 
de los rellenos manuales; la tecnologia 
apropiada para disponer los residuos sin incurrir 
en costos mayores de inversión y operación; y 
la ausencia de conocimiento acerca de cómo 
enfrentar el problema de la disposición final 
inadecuada de residuos. 

Residuo sólido: Es cualquier producto, 
materia o sustancia, resultante de la actividad 


humana o de la naturaleza, que ya no tiene más 


Capítulo 6 


función para la actividad que lo generó. 

Pueden clasificarse de acuerdo a su: origen (domiciliar, industrial, comercial, 
institucional, público), composición (materia orgánica, vidrio, metal, papel, plásticos, 
cenizas, polvos, inerte) y peligrosidad (tóxicos, reactivos, corrosivos, radioactivos, 
inflamables, infecciosos). 


21 OBJETIVOS 


Proveer de un sistema de relleno sanitario de manera a mitigar los impactos 
ambientales negativos. Que el RSM sea aceptado por la comunidad de Valenzuela y sea 
lo más económico posible para la misma, teniendo presente su sostenibilidad. Reducir el 
riesgo de incendios y accidentes debido a la acumulación de gases inflamables. Reducir 
los olores producidos por los residuos orgânicos. Reducir la contaminación debido a la 
descomposición de la materia orgánica 


3 | METODOLOGIA 


Se describe la zona del proyecto en los aspectos sujetos a impactos, y luego se 
elabora un Plan de Gestión Ambiental. 

Geomorfologia: el área es bastante homogénea en términos estructurales. El tipo 
de suelo es arcilloso. Suelo que es favorable para este tipo de proyecto desde el punto 
de vista de la impermeabilidad. Este suelo se utilizará como suelo de base del proyecto 
y además como relleno compactado diariamente sobre los residuos sólidos formando las 
denominadas celdas diarias. 

Relieve: el área está caracterizada por pendientes relativamente variables por lo que 
el drenaje superficial es bastante bueno y se observa escasa erosión. Siendo la pendiente 
del terreno natural 3%, presentando características favorables para la implementación del 
proyecto desde esta perspectiva. Clima: sub tropical húmedo, las temperaturas media 
mínima es de 13,8º Cy la máxima 33,2º C. Eco región: El proyecto esta ubicado en la eco 
región de Litoral Central. Descripción de las características del área de emplazamiento: 
La zona se caracteriza por la actividad agrícola ganadera. Cuenca Hidrográfica: la propiedad 
se encuentra ubicada entre los Aº Segovia Paso y Aº Paso Pe. Ubicados a 1.000 metros del 
lugar del emplazamiento del proyecto. - Tipo de vegetación (pastizal, arbustivo, arbóreo): 
Corresponde a la Eco región Litoral Central. Entre las especies arbóreas podemos citar 
Ybyrapyta, Curupay, mango. - Fauna: Se observan aves como, cardenal, tórtola, calandria 
para, hornero. - Asentamientos Humanos: El casco urbano de la ciudad está a una 
distancia de 3.000 metros. Inversión total: 120.000 U$S. Tecnologias y procesos que 
se aplicarán: La excavación se realizará en forma manual, el relleno sanitario también, 


diariamente contando con las herramientas y el personal apropiado para dichas prácticas. 
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El relleno sanitario es una técnica de disposición final de los residuos sólidos en el suelo 
que no causa molestia ni peligro para la salud o la seguridad pública; tampoco perjudica 
el ambiente durante su operación ni después de su clausura. Esta técnica utiliza principios 
de ingenieria para confinar la basura en un área lo más estrecha posible, cubriéêndola con 
capas de tierra diariamente y compactândola para reducir su volumen. Se contará además 
con el servicio de recolección domiciliaria de residuos con un recolector proveído por la 
misma municipalidad. Esta alternativa prevé los problemas que puedan causar los líquidos 
y gases producidos por efecto de la descomposición de la materia orgânica. Figuras 1 y 2: 





Figura 1: Lugar de emplazamiento del proyecto 





Figura 1: Toma de muestra de agua Ao Paso Pé 


4 RESULTADOS Y DISCUSIÓN 


Etapa 1: Disenio. En relación con la disposición final de RSM, de tres alternativas 
de tipos de rellenos sanitarios, a saber (se eligió la tercera debido al número de habitantes 
Po = 6.000, y a la tasa de crecimiento poblacional del 0,0% de la localidad y al volumen de 
basura diario generado 0,3kg / habitante dia): Pf=Po x (1+0,0)20=6.000 hab.(Proyección 
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población a 20 afios) Peso:365 dias x 6.000 hab. x 0,3 kg / habit dia =657.000 kg. / afo. 
Densidad relleno compactado 450 kg / m3; Volumen anual: 657.000 / 450 =1.460 m3/ afo; 
Volumen para 20 afos = 20x1.460 =29.200 m3 

El Relleno Sanitario Manual, en adelante RSM, es el elegido para esta comunidad, 
por su practicismo en su ejecución presupuestaria para el volumen y peso de basura diario 
de aproximadamente 4 mº / día y 1,8 ton / día, la tecnologia disponible en la zona y la 
población que cubrirá, siendo en este caso el menos oneroso para la misma. El terreno 
considerado para el RSM cuenta con un nivel freático de 6 metros de profundidad lo cual 
es apto para el proyecto que consistirá en: trincheras de 3 metros de profundidad con una 
base de 1,5 metros dejando un talud de pendiente 1: 2 para evitar el desmoronamiento del 
terreno adyacente y 1 metro entre trinchera y trinchera teniendo en cuenta las características 
del suelo y esta profundidad, se evitará la contaminación de las aguas subterrâneas y 
superficiales. 

El área de emplazamiento del proyecto requiere como mínimo: 

16.500 m2 para la zona de relleno sanitario propiamente dicho (29.200m83 / (3m 
x0,59) (coeficiente de aprovechamiento del terreno por el método de trinchera (ver corte de 
trinchera)= 16.497 m2) 

100 m2 para depósito, guardia, oficina del encargado y bafio del personal 

1.500 m2 para zona de vivero 

11.600 m2 para circulación de camiones y vehículos 

50.000 m2 de zona de amortiguamiento 

79.700 m2 de Área Total mínima del RSM 

-Precipitación pluvial en el área del Relleno Sanitario anual es de 1.219 mm 
(DGEEC-2006) - Escorrentia superficial (60%) - Evapotranspiración. - Humedad natural 
del RSM, 8 %.- Grado de compactación: relleno débilmente compactado K =0,35 (coef. 
de compactación) - Capacidad de campo (capacidad del suelo y del RSM para retener 
humedad: 40%). 

Coeficiente de permeabilidad del suelo: 


10-5 a 10-7 cm/s (muy baja infiltración); 10-7 cm/s (prácticamente impermeable) 
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DETALLE DEL AREA DE PROVECTO 
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Figura 3: Detalle del proyecto 


Se elaboró un Plan de Gestión Ambiental para el proyecto 


PLAN DE GESTION 
AMBIENTAL 


PROGRAMA DE 
MITIGACION DE 
IMPACTOS 


Proyecto de: 


I.Proyecto de manejo 


para evitar erosión (mov. 


amortiguamiento). 
3.Proyecto de control de 
vectores 

4.Proyecto de manejo de 
gases. 

5.Proyecto de manejo de 
líquidos lixiviados. 


PROGRAMA DE 
MONITOREO 


Proyectos de: 
1. Revisión de 
documentos. 
2. Evaluación de equipos 
operativos 
3. Revisión de limpieza y 
mantenimiento de 
equipos 
4. Monitoreo periódico de 
la calidad de aguas 
subterráneas. 
5. Monitoreo de gases 
6. Monitoreo de calidad 
de lixiviados luego de su 
tratamiento. 
7. Revisión de manuales 
de operación. 
8. Revisión de los equipos 
de seguridad. 
9. Revisión de los cursos 
de capacitación. 
10. Revisión de libros de 
quejas e incidentes 
ambientales. 





PROGRAMA DE 
SEGURIDAD Y 
MANTENIMIENTO 


Proyecto de: 
1. Equipos de protección 
individual (EPJ). 
2. Prevención de 
incendios y explosiones. 
3. Mantenimiento de 


equipos (preventivo y 
correctivo). 

4. Capacitación y 
entrenamiento. 

5. Mantenimiento de 
infraestructura. 

6. Sefialización. 


Diagrama: Plan de Gestión Ambiental 


Mitigación de impactos negativos del medio físico 


* Se tendrá en cuenta un plan para evitar erosión en los movimientos de suelos, 
que consistirá en colocar cobertura vegetal (pasto) en taludes. 


*— Arborización en la zona de amortiguamiento para servir de contención de polvo 
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y de eventuales particular en suspensión como ser pedazos de papel o plástico 
como además tendrá un efecto paisajístico. 


* | Dispositivos de drenaje de líquidos provenientes de la descomposición de resi- 
duos, que pasaran por una planta de tratamiento de este lixiviado, en una etapa 
primaria y otra secundaria hasta llegar a límites aceptables de contaminación. 
Esto además será acompariado de un monitoreo constante para identificar efi- 
ciencia del tratamiento. 


* -Quema controlada de gas metano proveniente de la descomposición de los 
residuos. 


Control de Impactos al hombre 


* | Plan de seguridad ocupacional, equipos de protección individual (EPI). Educa- 
ción y capacitación para los trabajadores y la sociedad para mejorar la disposi- 
ción de los residuos. 


* | Control de vectores, mediante insecticidas biodegradables. 


* | Plan de vacunación para los trabajadores. 


Jabla de Impactos y Medidas de Mitigación 


Etapa Impactos | lj Ea ita Era o | Responsable 


| Mingación | 


e O q a aan SÃO 


Construcción | | 
Movimiento de | Perdida de [Cobertura | | Supervisor 
suelo vegetaciôn | vegetal | Ambiental 
Construcción 

de ALLE SCS 


Excavación Accidentes | Plan de | Supervisor 


| 
| 


para laborales seguridad | Ambiental 
disposición 


final 





Tabla de Impactos y medidas de mitigación en la Construcción 


Coleção desafios das engenharias: Engenharia civil 4 Capítulo 6 


Tabla de Impactos y Medidas de Mitigación 


Etapa Impactos | Medidas de Responsable 
Mitigación 
EEE SE E 


Transporte, descarga y Generación de olores, | Localización basada en | Supervisor Ambiental 
compactación de residuos | polvo y vectores | el análisis de vento, 
barreras naturales 
Control de vectores 
Problemas de salud | Plan devacunacióny | Supervisor Ambiental 
ocupacional y | capacitación del 


poblacional personal 
Drenae de Liqudos Deterioro calidad del | Dispositivos de Supervisor Ambiental 


agua | encauzamento y 
tratamento 
Drenaje de Gases Descomposición de Dispositrros de Supervisor Ambiental 
resduos produce conducaión y quema 


metano potencialmente | controlada 
toxIcO F explossvo 





Tabla de impactos y mitigación durante la operación 


Tabla de Impactos y Medidas de 


4 Bis tes: Teto n 


= Impactos Medidas de Responsable 
Mitigación 


Clausura 


Arbornización Yv Aspecto Resntegrar Super VISOI 


cobertura paisajístico y efecto estético | Ambiental 


vegetal estético. Y funcional 
plantando 


Suelo 
arboles y 
cobertura 
vegetal( pasto), 
para evitar 


Eros. 





Tabla de impactos y mitigación etapa de clausura 


Desarrollar un plan de Emergencia 


* || Desarrollo de un plan de Respuesta a la Emergencia y entrenamiento de sus 
empleados en cómo usarlo. Ya que las emergencias son impredecibles, se 
debe preparar un Plan de Respuesta la Emergencia que refleje las condiciones 
de su entorno. 


* | Aldesarrollar dicho Plan de Respuesta a la Emergencia, considerar lo siguiente: 


* | Limite las acciones centralizando las actividades alrededor de la Emergencia. 
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* | El plan debe basarse en un número mínimo de empleados presentes en el 
relleno sanitario. 


* | Elplan debe estar expuesto y claramente visible para conocimiento de todos. 


* | El entrenamiento de su personal en la ejecución del plan le asegura un alto 
grado de éxito en el manejo de emergencias, de manera que entrene, entrene 
Y siga entrenando a su personal. 


Obligaciones del supervisor ambiental 


* | Promover constantemente la seguridad y la higiene. 
* | Control periódico y monitoreo de los drenajes líquidos y gases. 
* | Inventario diario de la cantidad de residuos que Ilegan. 


* | Capacitación del personal cada 12 meses, sobre riesgos existentes de acuerdo 
a funciones que desempefa. 


* Vigilancia a la salud según establece IPS. 

* | Manual de equipos. 

* | Programa de revisión de equipos e instalaciones. 
* | Simulacro de accidentes 

* | PAE (Plan de Atención a Emergencias) 


* | Plan de respuesta ante accidentes de personal. 


Obligaciones de los trabajadores 


* Estar informado de los riesgos existentes en la planta. 

* | Participar en la capacitación proporcionada. 

* | Acceder a la práctica de exámenes médicos. 

* | Desarrollar las actividades de acuerdo al manual de operaciones. 


* Utilizar los equipos de protección individual de acuerdo al trabajo a ejecutar, 
mantenerlos en condiciones de uso y limpieza e informar al Capataz cuando 
debe reponerse. 


* Seguir las Instrucciones del Supervisor Ambiental o Capataz, para realizar ac- 
tividades peligrosas. 


* | Participar en las capacitación de Seguridad e Higiene. 
Proyecto de capacitación 
Se debe instrumentar bajo el concepto de: 


. Decir 
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* Mostrar 
* Indicar 
* - Corregir 
*  Supervisar 
Plan de Monitoreo 
El objetivo del plan es establecer para cada fase del proyecto o actividad: 
* || Los componentes del medio ambiente que serâán objeto de medición y control 
* | Elimpacto ambiental asociado 
* | La ubicación de los puntos de control 


* | Los parâmetros que serán utilizados para caracterizar el estado y evolución de 
dicho componente 


* | Los niveles cuantitativos o limites permitidos o comprometidos 
* La duración y frecuencia del plan de monitoreo para cada parâmetro 
* | El método o procedimiento de medición de cada parâmetro 


* | El plazo y frecuencia de entrega de los informes del plan de monitoreo a los 
organismos competentes. 


* | Laindicación del organismo competente que recibirá dicha documentación 
* | Cualquier otro aspecto relevante 
Los componentes a ser controlados: 
* | Calidad del aire 
* | Calidad de las aguas subterrâneas 
* | Calidad de las aguas superficiales 
* | Control del asentamiento del relleno, agrietamiento y erosión de superficies 
* | Control de los sistemas de drenaje (lixiviados y aguas de Iluvia) 
«| Calidad del efluente (lixiviados) 
Frecuencia 


* | El muestreo y análisis de la calidad de las aguas superficiales se efectuara con 
la siguiente frecuencia: 


* Antes de iniciar la operación: una vez, análisis de todos los parâmetros. 
* | Durante la operación: semestral 
* | Después del cierre: una vez, análisis de los parâmetros de contingencia 


. Durante el abandono: anualmente, análisis de los todos los parâmetros durante 
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tres anos. 


Parámetros de control 


- PH 
* | Temperatura 
* Aceitesy Grasas 
*- | DQO (Demanda Química de Oxigeno) 
* Manganeso 
*- Conductividad 
*  Nitrógeno Total 
*  Cloruros 
* Sulfatos 
Monitoreo del lixiviado 


* | Los líquidos lixiviados generados por el RSM que tienen su origen en dos fuen- 
tes, fundamentalmente el agua de Iluvia infiltrada y las internas, liquido propio 
de los residuos y el generado en los procesos de descomposición de la materia 
orgânica presente. Ambos pasan a través de los residuos depositados, diluyen- 
do y arrastrando diversos componentes de la basura. El contacto entre la fase 
liquida enunciada y la sólida compuesta por basura, produce un líquido final que 
contiene una mezcla de elementos orgánicos e inorgánicos en diferentes for- 
mas según la movilidad de cada uno (solución, emulsión, sólidos arrastrados), 
que deben ser analizados. 


Control de Asentamiento del relleno, Agrietamiento y Erosión de Superficies 


* | Serân controlados en forma visual y a través de evaluaciones topográficas alti- 
métricas cada afio. 


Sitios de Monitoreo 


* - Superficie del relleno sanitario. 


Frecuencia 


* ' Mensualmente durante la operación y anualmente durante tres aos después 
del cierre. 


Parámetros 


* Se realizara el monitoreo considerando alturas, pendientes, posibles impactos 
por erosión y formación de grietas. 
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Control de los Sistemas de Drenaje ( Lixiviados y Aguas de Lluvia) 


Metodologia 


* | Las tuberías de drenaje cuentan con câmaras de inspección y control, en forma 
visual. 


* | Las canaletas de agua de Iluvia serán inspeccionadas en forma visual. 


Sitios de Monitoreo 


*  Tuberías de drenaje para conducción de líquidos percolados y canaletas de 
drenaje de Iluvia. 


Frecuencia 


* | Lafrecuencia de inspección será: 
* | Durante la operación : Semanal (control visual) 


* | Durante el abandono: Mensual y después de Iluvias intensas. 


Parámetros 


* | Se monitoreara el estado del sistema de drenaje, inspeccionando tuberías y 
canaletas en busca de taponamientos y roturas. 


5 | CONCLUSIONES 


Es importante que en se tenga a la ciudadanía concientizada sobre la limpieza de la 
comunidad y la importancia de los proyectos asociados de recolección de basura, limpieza 
de calles y lugares públicos, el reciclaje y la utilización del RSM Piha Porá para tener 
una comunidad saludable y orgullo p/habitantes. Se accede a la propiedad, por la ruta 
empedrada Itacurubi-Valenzuela y por un camino vecinal a ser construido. 

Una ventaja del relleno sanitario como método de tratamiento de residuos es la 
posibilidad de recuperación de áreas ambientalmente degradados por explotación de 
canteras, así como terrenos considerados improductivos o marginales, que se pueden 
convertir en parques y plazas para el usufructo de la comunidad una vez cumplido con los 
trabajos de clausura. 

El relleno sanitario manual es de baja inversión y mantenimiento, para lo que 
además debe ser controlado, hacer cumplir los diferentes programas de mitigación de 
impactos, cumplir con los programas de monitoreo para servir de base indicativa de la 
eficiencia de los programas de mitigación, contar con un monitoreo de línea de base que 
sirva para identificar alteraciones en las aguas subterrânea y aguas superficiales próximas, 


para identificar posibles alteraciones del entorno. 
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